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This  study  has  been  conducted  by  the  Pomona  Operation  of  General 
Dynamics  Corporation  for  the  Naval  Air  Systems  Command  under  Contract 
N00019-69-C-0319. 

The  principal  investigator  for  the  study  has  been  Mr.  John  G. 
Maloney.  The  experimental  work  and  much  of  the  analytical  work  was 
carried  out  by  Mr.  Michael  T.  Shelton.  Dr.  George  Lasker  served  as  a 
staff  consultant  and  provided  the  finite  element  analysis  for  the 
correlation  analysis.  Mr.  Kenneth  L.  McIntyre  developed  the  joint 
compliance  extraction  technique  and  the  joint  spring  coupling  analysis. 
The  direct  technical  supervisor  has  been  Mr.  David  A.  Underhill, 
Structural  Dynamics  Section  Head. 

The  Naval  Air  Systems  Command  technical  monitor  has  been 
Mr.  George  P.  Maggos. 
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Section  1.0 
SUMMARY 


The  initial  results  of  a  study  of  the  structural  dynamic 
properties  of  tactical  missile  joints  are  presented.  The  scope  of  this 
phase  of  the  investigation  has  included: 

1.  A  review  of  the  types  of  mechanical  joints  in  common  usage 
and  the  methods  employed  in  estimating  and  representing 
joint  load/deflection  characteristics  in  structural  dynamic 
response  studies. 

7.  A  parametric  study  to  identify  and  illustrate  the  controlling 
mechanisms  and  relative  importance  of  joint  and  airframe 
structural  properties  and  geometries. 

3.  The  preliminary  evaluation  of  a  method  for  extracting  joint 
compliances  from  measured  modal  data. 

4.  A  simplified  test  series  using  idealised  models  to  explore 
and  illustrate  the  effects  of  load  path  discontinuities  and 
dissymmetries . 

Based  on  comments  and  replies  received  in  an  industry  survey 
(included  as  an  appendix  to  this  report),  tactical  missile  joints  ere 
generally  represented  in  dynamic  analysis  by  equivalent  rotational  springs 
selected  by  trial  and  error  to  match  measured  response  characteristics. 
Most  respondents  cited  significant  reductions  in  airframe  flexural  mode 
frequencies  due  to  joint  deflections,  indicating  substantial  losses  in 
local  airframe  stiffness.  Using  a  classification  basis  proposed  in  a 
NASA  study,  a  joint  considered  to  be  "moderate"  in  compliance,  for 
example,  is  shown  to  result  typically  in  a  ninety-five  percent  loss  in 
effective  airframe  stiffness  for  a  reference  span  of  half  a  body  diameter. 

Parametric  analysis  and  test  results  illustrate  the  effects  of 
joint  location  and  rating  on  airframe  dynamic  response  characteristics 
and  show  the  sensitivity  of  joint  compliance  to  load  path  dissyracetry. 

Load  path  discontinuities  and  dissynnetries  are  concluded  to  be  a  major 
contributor  to  joint  compliance. 

A  method  of  steepest  ascent  applied  to  the  problem  of  extracting 
joint  compliance  from  measured  modal  data  is  developed  and  shown  to  have 
considerable  promise.  Direct  analysis  of  joint  characteristics  using 
finite  element  modeling  techniques,  applied  in  a  test  case  for  correla¬ 
tion  purposes,  similarly  shows  significant  potential  as  an  analytical 
tool.  The  status  of  the  study  at  the  completion  of  the  first  phase  is 
reviewed  and  the  objectives  of  the  second  phase  outlined. 
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Section  2.0 
INTRODUCTION 


The  trend  with  high  performance  missiles  as  with  aircraft  is 
toward  structural  design  rsquirements  in  which  airframe  stiffness  plays 
an  increasingly  important  role.  The  source  of  airframe  stiffness 
requirements  can  stem  from  a  variety  of  system  design  considerations 
including  static  and  dynamic  aeroelastlc  stability  margins,  airframe 
aeroelastic  coupling  with  guidance  and  control  systems,  and  structural 
dynamic  loads  and  response  induced  by  logistic  and  flight  environments. 

Experience  has  shown  that  many  of  the  mechanical  joints  conmtonly 
employed  in  tactical  missiles  to  serve  modular  design  requirements  can 
result  in  substantial  and  often  unpredictable  reductions  in  the  stiff¬ 
ness  of  the  primary  structure.  In  the  absence  of  reliable  analysis 
methods  for  estimating  joint  effects  on  airframe  stiffness,  conn  on 
practice  is  to  rely  on  experimental  data  for  definition  of  joint  proper¬ 
ties.  The  shortcoming  of  this  approach,  however,  is  that  data  obtained 
for  a  particular  joint  design  on  a  given  missile  often  cannot  be  extra¬ 
polated  with  any  confidence  to  a  different  airframe  design  or  even,  in 
many  cases,  to  a  different  location  on  the  3ane  airframe.  The  lack  of 
reliable  methods  for  predicting  the  load/deflection  behavior  of  mechani¬ 
cal  joints  is  a  limitation  of  increasing  concern  in  the  development  of 
efficient  structures  for  advanced  tactical  missile  airframes  and  in  the 
early  assurance  of  structural  integrity  in  the  design  development  phase. 

The  study  described  in  this  report  represents  the  first  phase  in 
a  basic  investigation  of  the  structural  dynamic  properties  of  tactical 
missile  joints  being  undertaken  by  the  Pomona  Division  of  General 
Dynamics  for  the  Naval  Air  Systems  Command. 

One  of  the  initial  tasks  has  been  a  compilation  of  the  types  of 
mechanical  joints  in  common  usage  together  with  a  review  of  the  current 
methods  employed  in  estimating  and  representing  joint  load/deflection 
characteristics  in  structural  dynamic  response  s.'tudies.  This  effort  has 
included  a  literature  search  as  well  as  an  aerospace  industry  survey  and 
is  summarised  in  Sections  3,0  and  4.0. 

Following  this,  a  parametric  study  was  undertaken  to  investigate 
the  controlling  mechanisms  and  relative  importance  of  joint  structural 
properties  by  examining  their  effects  on  the  structural  dynamic  response 
of  representative  missile  structures.  The  results  of  this  study  are 
presented  and  discussed  in  Section  5.0  with  the  first  part,  Section  5.1, 
devoted  to  the  effects  of  variations  in  joint  compliance,  location  and 
number,  and  the  second  part.  Section  5.2  concerned  with  joint  dlssymetry 
and  elastic  coupling  behavior, 
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Section  6.0  presents  a  preliminary  evaluation  of  an  analysis 
procedure  for  extracting  joint  compliances  from  measured  modal  response 
data.  The  assumption  is  made  that  in  most  cases  the  basic  airframe 
stiffness  and  weight  distributions  are  adequately  defined,  with  the 
prLmary  unknown  being  the  airframe  joint  compliances.  Given  a  set  of 
measured  mode  shapes  and  frequencies,  the  objective  is  to  devise  a 
rapid  means  for  converging  on  the  effective  joint  compliances.  The 
accuracy  and  limitations  in  the  current  form  of  the  analysis  are  reviewed. 

A  limited  test  program  described  in  Section  7.0  has  been  accom¬ 
plished  using  highly  simplified  models  to  illustrate  and  verify  some  of 
the  parametric  analysis  results.  The  models  employ  section  property 
discontinuities  and  dissymmetries  Intended  to  simulate  mechanical  Joint 
behavior. 

Section  8.0  presents  a  status  review  for  the  first  phase  of  this 
investigation  and  outlines  the  scope  of  the  study  planned  for  the  second 
phase. 
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Section  3.0 

THE  STRUCTURAL  DYNAMIC  ROLE  OF  TACTICAL  MISSILE  JOINTS 


Mechanical  joints  are  a  common  feature  of  tactical  missile  struc¬ 
tures  used  to  join  major  airframe  sections  such  as  guidance,  control, 
ordnance,  autopilot,  and  propulsion.  The  great  variety  in  the  types 
of  joints  used  for  these  purposes  is  illustrated  in  the  results  of  an 
aerospace  industry  survey  presented  in  the  Appendix  to  this  report. 
Sketches  of  some  typical  Joints  in  comnon  usage  in  tactical  missile 
airframe  design  are  shown  in  Figure  3-1. 

Although  both  separable  and  non-saparable  joints  are  included  in 
the  survey,  the  ones  of  primary  interest  to  this  study  are  of  the 
separable  design.  Non-aeparable  joints  designed  for  permanent  attach¬ 
ment  of  sections  tend  to  be  stiffer,  stronger,  consistent  and  more 
predictable  in  their  properties  by  virtue  of  being  welded,  bonded  and/ 
or  riveted  with  many  fasteners.  Separable  joints,  however,  especially 
those  designed  for  ease  of  assembly  and  disassembly  under  field  and 
depot  conditions,  exhibit  behavior  under  dynamic  loading  which  is  often 
difficult  to  predict  and  sensitive  to  many  parameters  which  cannot 
readily  be  specified  or  controlled.  Torque  values,  for  example,  are 
frequently  specified  for  joint  fasteners  in  an  attempt  to  control  inter¬ 
face  preloads.  Measurements  made  with  Instrumented  bolts,  hcwsver, 
have  shown  that  the  axial  load  In  bolts  torqued  to  the  same  value  can 
vary  significantly  due  to  dimensional,  frictional,  and  thermal  variations. 

The  characteristics  of  airframe  joints  that  Influence  the  structural 
dynamic  behavior  of  tactical  missiles  can  be  categorised  as  follows: 

1.  Compliance 

2.  Energy  dissipation 

3.  Energy  transfer 

3.1  COMPLIANCE 

The  most  pronounced  characteristic  of  the  typical  airframe  joint  in 
a  tactical  missile  is  ita  local  compliance.  A  joint  constitutes  a 
disturbance  in  load  p3th  which  can  result  in  substantial  losses  in 
effective  stiffness  in  the  vicinity  of  the  joint.  The  consequences  of 
joint  stiffness  losses  3re  well  illustrated  by  examining  the  bending  modes 
for  an  actual  tactical  missile.  The  missile  used  in  this  example  is  the 
Standard  ARM.  The  type,  location,  and  estimated  compliance  of  the  six 
principal  airframe  joints  are  shown  in  Figure  3-2.  The  bending  modes  for 
the  airframe  were  computed  using  a  conventional  lumped  parameter  beam 
representation  in  a  modified  Holaer^Iyklestad  modal  analysis  method 
(Reference  7).  Values  for  the  joint  compliances  were  derived  by  a  trial 
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and  error  matching  of  the  airframe  modes  measured  in  ground  vibration 
testing.  The  first  three  bending  mode  frequencies  for  this  airframe 
with  and  without  joint  compliance  are  tabulated  below  and  the  corres¬ 
ponding  mode  shapes  shown  in  Figure  3-3.  In  this  example,  the  influence 
of  internal  appendages  and  the  control  surfaces  is  purposely  omitted  in 
order  to  show  the  direct  effect  of  the  airframe  joint  compliance. 


Table  3-1 


JOINT  EFFECTS  ON  STANDARD  ARM  MODES 


Mode 

No. 

Mode  Frequency  (Hz) 

Freq. 

Loss 

% 

Stiffness 

Loss 

1 

Without 

Joints 

With 

Joints 

1 

62 

51 

18 

32 

2 

159 

115 

28 

48 

3 

302 

206 

_ 

32 

53 

The  estimated  equivalent  loss  in  airframe  generalized  stiffness  for 
each  mode  is  based  on  the  approximate  assumption  of  no  change  in  general¬ 
ized  mass: 

/JLY  *  .  **.  « 

K  *1  /  Kj 


This  assumption  is  not  strictly  speaking  true  since  changes  in  mode 
shape  are  obvious,  but  it  does  give  a  fair  indication  of  the  powerful 
influence  that  the  joints  have  on  the  airframe  stiffness. 

Table  3-2  presents  for  comparison  purposes  similar  information 
obtained  in  the  Industry  Survey  (Appsndlx  I)  tor  several  other  tactical 
missiles.  In  this  case,  data  only  on  the  first  mode  are  shown.  The 
equivalent  stiffness  loss  for  the  airframe  has  again  been  estimated  on 
the  same  basis  as  previously  discussed. 

The  significant  point  is  that  substantial  joint  compliance  effects 
are  evident  on  many  representative  tactical  missile  airframes.  In  view 
of  this,  it  is  clear  that  accurate  estimates  of  joint  compliance  are 
likely  to  be  of  critical  importance  in  predicting  airframe  response 
characteristics . 
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Table  3-2 

ESTIMATED  JOINT  EFFECTS  -  MISSILE  FIRST  MODE  CHARACTERISTICS 


Missile 

Number  of 
Joints 

1st  Mode 
Freq.  Lose 

% 

Stiffness 
Los  8 

% 

Sidewinder 

A 

7 

13 

SRAM 

6 

17 

31 

Standard  (MR) 

7 

17 

31 

Standard  (ER) 

7 

18 

33 

Phoenix 

10 

29 

49 

3.2  ENERGY  DISSIPATION 

Tactical  missile  joints  clearly  play  a  major  role  in  dissipating 
vibratory  energy.  This  conclusion  is  unavoidable  considering  the  nature 
of  tactical  missile  structures  and  the  fact  that  a  substantial  portion 
of  the  total  potential  energy  in  the  lower  modes  of  interest  is  stored  in 
the  joint  load/deflection  characteristics.  Missile  airframe  structure 
between  joints  typically  consists  of  a  simple  single  element  cylindrical 
section  (as  opposed  to  built  up  structure)  in  which  the  only  source  of 
damping  is  that  inherent  in  the  material  properties.  Connections  associated 
with  internal  appendage  mountings  of  course  constitute  another  source  of 
energy  dissipation,  but  are  not  believed  to  be  a  major  contributor  to  air¬ 
frame  damping  in  the  lower  modes. 

The  precise  nature  of  energy  dissipation  in  mechanical  joints  is  not 
well  understood  but  is  believed  to  involve  two  basic  mechanisms: 

1.  Sliding  friction 

2.  Gas  pumping 

Sliding  friction  is  suspected  to  be  the  dominant  source  in  the  low  to 
moderate  frequency  range  with  gas  pumping  effects  becoming  important  at 
higher  frequencies  (Ref.  6). 

It  is  interesting  to  note  that  the  modal  damping  of  missile  struc¬ 
tures  is  generally  in  the  same  domain  as  aircraft  structures  (at  least 
in  the  lower  modes,  typically  ranging  from  1  to  3  percent  of  critical) 
although  their  structural  configurations  are  considerably  different. 
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Aircraft  structural  damping,  however,  tends  to  be  distributed  more 
uniformly  over  che  structure  rather  than  localized  at  a  few  major  joints 
as  in  the  case  of  missile  structures.  Missile  structural  dynamic 
response  in  consequence  often  exhibits  abrupt  and  sizable  phase  shifts 
across  joints. 

3.3  ENERGY  TRANSFER 

Missile  airframe  joints,  in  addition  to  their  compliance  and 
damping  characteristics,  are  responsible  for  two  important  forms  of 
energy  transfer: 

1.  Elastic  coupling  of  missile  response  coordinates. 

2.  Energy  transfer  in  the  frequency  domain  through  joint  interface 
non-linearities  and  impact. 

The  elastic  coupling  behavior  of  joints  is  believed  to  result  from 
a  non-axially  symmetric  distribution  of  load  paths  through  a  joint 
interface.  A  conceptual  model  of  this  is  developed  and  applied  in 
Section  5.2  of  this  report.  A  representative  example  of  cross  plane 
response  attributed  to  joint  elastic  coupling  on  an  actual  tactical 
missile  is  presented  in  Figure  3-4.  The  data  describes  the  angular 
response  in  pitch  (in-plane)  and  yaw  (cross-plane)  at  the  autopilot  rate 
gyro  station  for  a  unit  force  at  the  control  surface  station  driving  the 
pitch  plane.  In  this  example,  the  cross -plane  response  is  within  2  dB 
of  the  in-plane  response  in  both  of  the  first  two  dominant  modes. 

The  transfer  of  energy  in  the  frequency  domain  attributable  to  air¬ 
frame  joints,  is  a  result  of  joint  stiffness  non-linearities.  This 
characteristic  manifests  itself  in  two  forms:  as  harmonic  distortion  of 
the  response  with  certaia  forms  of  stiffness  non-linearities, -and  as  a 
high  frequency  noise  source  when  free  play  results  in  impact  of  the  con¬ 
tacting  surfaces  of  the  joint. 

An  example  of  the  harmonic  distortion  response  which  an  airframe 
joint  can  produce  is  taken  from  Reference  5,  These  results  were  obtained 
during  an  investigation  of  the  effect  of  out -of -tolerance  thread  lead  on 
the  coupling  nut  used  in  the  airframe  joint  shown  in  Figure  1-15.  The 
tests  were  conducted  using  a  dummy  forward  section  attached  to  a  vibra¬ 
tion  fixture  via  the  co  >ling  nut  joint.  The  assembly  was  excited  with 
4g  peak  sinusoidal  vibration  in  a  lateral  direction  at  30  Hz,  which  was 
approximately  two  octaves  in  frequency  below  the  first  lateral  resonance 
of  the  assembly.  It  was  found  that  the  presence  of  a  set  screw  in  the 
coupling  nut  (used  to  lock  the  assembly)  had  a  profound  influence  on  the 
assembly  response  characteristics.  The  lateral  response  acceleration 
time  history  was  monitored  at  a  joint  in  the  dummy  guidance  section  well 
removed  fran  the  joint.  The  acceleration  traces  are  presented  in  Figure 
3-5a  for  one  of  the  coupling  nut  specimens  with  the  set  screw  removed  and 
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with  the  set  3crew  properly  installed.  T.t  can  be  seen  that  the  accelera¬ 
tion  response  is  a  very  clean  sinusoid  when  the  set  screw  is  properly 
installed  but  that  harmonic  distortion  of  the  response  occurs  when  the 
set  screw  is  not  installed. 

An  example  of  the  generation  of  high  frequency  noise  due  apparently 
to  free  play  in  a  joint  can  be  seen  in  Figure  The  test  condi¬ 

tions  are  exactly  the  same  as'  those  described  above.  The  coupling  nut 
has  been  replaced  with  a  second  sample.  The  absence  of  the  set  screw 
results  in  an  acceleration  response  with  considerable  high  frequency 
energy  content.  In  this  instance  the  harmonic  distortion  and  high 
frequency  noise  occurred  only  when  the  joints  were  improperly  installed, 
that  is  when  the  set  screws  ware  not  in  place.  These  examples  have  been 
cited  only  to  show  In  a  qualitative  sense  the  phenomena  of  airframe 
joint  energy  transfer  characteristics  in  the  frequency  domain. 


8 


POMONA  OPERATION 
(01  E-M* 


ggjjjg*k  BVWAMrg 
Electro  Dynamic  Division 


re'/rSHJ&JE  3-/ 

rv/o/c^d.  At /&■&/*.£  so/a/7~s 


WrEd.dSrED 


3OAS0JZO 


<e*sje7Z£o 


s sosast^ 


^e««a£Z) 


>f^e*rA(/ 


dA?ASO 


3*S£> 94* 


JT4*/4^4!ac  Sosas srs" 


POMONA  Of  {RATION 

;o;  i-hs 


9 


aWNgRAL  DYNAMICS 

Eiactro  Dynamic  Division 


3~Z 

&7V94A 0**0  **M  >/<3/X/T  0*7* 


Z*tS£L  £^£9USOL 
£*/*>z*so/x  ) 

AS/AVG  Z-/9 

V-ccov^"  jr-zo 

CLa*A4/0 

T-Z / 

TJEA/540A/  JT-ZZ 


SMB**  z-ZS 

72EVS/OA/  Z~ZZ 


aoA*Au*A/cje 


M/-/.A 


0.70 


iT.O 


Z.O 


O./ 


/.O 


10 


POMONA  OPERATION 
i  O)  6-S«S 


I'wrramiCT’WwniiP 


agNBF 

ESvCifu  L 


I 

I* 


% 


POMONA  0 


fNBRAL  DYNAMICS 

Electro  Dynamic  Division 


3-+ 

<3A=-  <£*0SS  /O^ASJE  *ejers/*eDAS'SJe: 


0 


/Ay  A^CAPA/^e  AZJZSA^&AS&je 
&eass  A?A./9A/r£'  /eJS.SiBOA/SJE' 


3 -Sol 

s£X./9M/4c<€'  <3^  /v&eAt£>AS/c  &'irr<2&r/d>AS 
GcJve:  7&  jo/A/r 


W/r//  xxr-  zsrcxTje^/ 


W'movr-  ST/-  £-&ej^u/ 


ats&cucjz :  3 -^"6 

s£7L/9A** 5CJC  ^  AWST  <S^GO€»?7515A' 
z>44«r  ro  sr/je^/estAtc  joftsr- 


\ 

% 


W/7W  s^r  src&ev/ 


W/r/W7-  s-^y-  cs&PjEW 


/3 


qiBNBWAL  DYNAMIC* 

Electro  Dynamic  Division 


Section  4.0 

STRUCTURAL  DYNAMIC  MODELING  OF  AIRFRAME  JOINTS 


4.1  AIRFRAME  REPRESENTATION 

A  common  tpproach  to  modeling  the  structural  dynamic  characteristics 
of  tactical  missiles  consists  of  developing  a  non-uniform  beam  model 
using  a  lumped  parameter  representation  for  the  distributed  stiffness 
and  weight.  Mode  shapes  and  frequencies  are  computed  for  this  representa¬ 
tion  by  a  variety  of  standard  modal  analysis  procedures  in  order  to 
provide  generalized  coordinates  for  use  in  all  types  of  structural 
dynamic  and  aeroelastic  response  studies.  The  adequacy  of  the  non-uniform 
beam  model  is  of  course  dependent  upon  structural  details  peculiar  to  the 
individual  missile  design,  with  particular  emphasis  in  the  present  study 
focused  on  the  load/deflection  behavior  of  the  missile  joints. 

The  development  of  the  lumped  parameter  representation  for  the  air¬ 
frame  weight  and  stiffness  distribution  (with  the  exception  of  the  joint 
compliance  characteristics)  is  usually  a  straightforward  process  involv¬ 
ing  the  simple  geometric  and  material  properties  of  the  structure. 
Mechanical  joints,  however,  constitute  a  more  difficult  and  less  straight¬ 
forward  modeling  task. 

The  results  of  the  industry  survey  indicate  that  load/deflection 
characteristacs  of  mechanical  joints  in  missile  airframes  are  most 
frequently  represented  by  a  flexural  or  rotational  spring.  Some  beam 
modal  analysis  programs  directly  admit  the  assumed  joint  stiffness  or 
its  reciprocal,  the  joint  compliance,  while  other  programs  require  that 
the  joint  effects  be  accounted  for  by  reducing  the  airframe  stiffness  In 
the  local  region  of  the  joint.  One  advantage  of  this  latter  procedure 
from  a  conceptual  standpoint  is  that  the  amount  of  local  stiffness 
reduction  (required  to  account  for  the  joint  compliance  contribution) 
provides  a  direct  basis  for  judging  the  significance  of  the  joint  rela¬ 
tive  to  the  basic  airframe  stiffness.  A  simple  relationship  between 
these  two  common  methods  for  representing  joint  compliance  can  be  shown 
by  equating  the  net  change  in  slope  per  unit  moment  over  the  reference 
airframe  length  selected  for  the  local  stiffness  reduction  as  follows: 


Reduced  Stiffness  Rotational  Spring  Model 
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where: 


/ 

&Ai 


^*4» 


reference  length  over  which  airframe  stiffness  is 
assumed  to  be  reduced. 

total  bending  slope  change  per  unit  moment  over  the 
reference  length. 

basic  (unmodified)  airframe  section  properties  in 
the  region  of  the  joint. 

effective  (reduced)  airframe  section  properties  in 
the  region  of  the  joint. 

material  flexural  modulus  of  elasticity. 

rotational  compliance  attributable  to  the  joint 


Combining  these  expressions  and  letting  JCg  denote  the  stiffness 
loss  or  reduction  ratio, 


In  comparing  the  two  methods  for  joint  representation,  the  selec¬ 
tion  of  the  effective  length  for  the  joint  influence  cn  airframe  stiff¬ 
ness  Is  somewhat  arbitrary.  Both  models  are  idealizations  which  over¬ 
simplify  the  structural  deformation  in  the  local  region  of  the  joint. 
Additionally,  the  deflection  models  for  the  two  methods  are  equivalent 
only  in  the  sense  of  matching  slopes  for  applied  moments.  The  signifi¬ 
cance  of  this  difference  is  usually  small  if  the  effective  length  is 
on  the  order  of  the  body  diameter  or  less.  Modal  analyses  were  conducted 
for  the  tactical  missile  example  given  in  Section  3.1  using  both  methods 
of  joint  representation  for  comparison.  An  effective  length  of  half  a 
body  diameter  for  joint  induced  stiffness  loss  was  used  in  this  compari¬ 
son,  and  virtually  no  differences  in  either  frequencies  or  mode  shapes 
were  noted  for  the  first  three  modes.  The  choice  of  half  a  body 
diameter  for  a  reference  length  in  judging  joint  effects  appears 
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reasonable  and  basing  the  reference  length  on  body  diameter  will  be 
shown  to  be  useful  from  a  dimensional  analysis  and  scaling  standpoint 
in  parametric  studies. 

4.2  DETERMINATION  OF  JOINT  COMPLIANCE 


T!  j  more  difficult  question  in  modeling  tactical  missile  airframe 
joints  is  not  how  but  what  to  model.  Compliance  estimates  based  on  local 
section  properties  and  fastener  elastic  characteristics  are  invariably 
uncon3ervative  (low).  If  representative  missile  hardware  exists,  one 
approach  is  to  perform  either  static  of  dynamic  load/deflection  tests  in 
which  joint  compliances  are  directly  measured.  Another  common  experimen¬ 
tally  based  technique  is  to  match  by  trial  and  error  a  measured  set  of 
mode  shapes  and  frequencies,  assuming  joint  compliances  to  be  the  only 
unknowns  in  the  modal  analysis  representation.  An  attempt  to  automate 
this  approach  is  described  in  Section  6.0. 

In  some  instances  a  ver  similar  joint  has  been  used  in  ar»  earlier 
application  and  an  extrapolation  of  its  compliance  characteristics  can 
be  made  with  acceptable  confidence.  This  is  not  always  th^  case,  how¬ 
ever,  and  there  are  many  examples  where  a  given  joint  design  behaves 
quite  differently  on  different  airframes  or  in  different  locations. 

Lacking  experimental  data  or  the  opportunity  to  test  representative 
hardware,  which  is  usually  the  case  in  pre-design  or  early  design  studies, 
the  analyst  is  confronted  with  the  task  of  making  a  best  judgement 
estimate  of  joint  compliance  characteristics.  One  approach  to  joint 
compliance  estimation  being  used  by  several  of  the  respondents  to  the 
industry  survey  is  that  outlined  by  Alley  and  Leadbetter  of  NASA  Langley 
in  Reference  1.  In  this  reference  Alley  and  Leadbetter  present  order  of 
magnitude  relationships  between  joint  flexural  compliance  and  joint  type 
which  they  derived  from  launch  vehicle  test  data.  The  compliance 
characteristics  are  classified  from  excellent  (smail  compliance)  to  loose 
(large  compliance)  according  to  the  table  shown  in  Figure  4-1  which 
covers  the  various  types  of  joints  examined  in  their  study. 


The  relationship  of  the  magnitude  of  joint  compliance  to  the  joint 
classification  is  shown  in  the  curves  in  Figure  4-1.  Alley  and  Lead- 
better  established  these  curves  based  on  10  measured  values  of  compliance 
and  the  assumption  that  the  compliance  is  inversely  proportional  to  the 
third  power  of  the  airframe  diameter.  The  curves  given  in  their  figure 
may  be  approximated  by  the  following  relationship. 


where 


C  =  flexural  compliance,  radians /inch -pound 

D  =  diameter  of  the  airframe  at  the  joint 
location,  inches 
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A  =  compliance  coefficient,  radians/inch- 
pound 

The  magnitude  of  the  compliance  coefficient  A  is  a  function  of  compliance 
classification  and  is  found  in  the  following  table. 


COMPLIANCE 

COMPLIANCE  COEFFICIENT  -  A 
(Radians/Inch-Pound) 

CLASSIFICATION 

NOMINAL 

RANGE 

Excellent 

o 

♦— t 

1 

/-v 

o 

r-4 

r-f 

<  3(10) “l0 

Good 

1(10)‘9 

3(10) _1G  to  3(10)“ 

Moderate 

1(1G)“8 

3(10) to  3(10)' 

Loose 

1(10) 

>  3(10)'® 

4.3  JOINT  COMPLIANCE  COMPARISONS 

The  results  of  the  industry  survey  provide  a  basis  for  comparing 
some  of  the  compliance  values  which  have  been  assigned  to  various  tacti¬ 
cal  missile  joints.  Using  the  Alley/Leadbetter  proposed  classification 
criteria  as  a  framework  for  comparison.  Figure  4-2  presents  joint  com¬ 
pliance  estimates  for  four  tactical  missiles  ranging  from  2.75"  to  18"  in 
diameter.  Although  only  data  for  single  joints  are  shown- for  Sidewinder 
and  Thor,  the  range  of  Standard  ARM  joint  compliances  is  seen  to  vary 
nearly  two  orders  of  magnitude.  An  additional  comparison  with  the  "NASA" 
rating  criteria  using  the  Standard  ARM  as  an  example  is  given  in  Figure 
4-3.  In  this  comparison,  an  additional  parameter,  stiffness  reduction 
ratio  (A*)  equivalent  to  the  joint  compliance  has  been  computed  using 
the  expression  for  developed  in  Section  4.1.  It  is  particu.arly 
revealing  to  note  that  for  this  airframe,  even  the  joints  considered 
"good"  represent  a  local  (half  body  diameter)  stiffness  loss  of  nearly 
70  per  cent,  while  the  "moderate"  joints  approach  95  percent  stiffness 
loss. 


The  applicability  of  this  conclusion  to  other  missiles  is  shown  in 
Figure  4-4  which  is  based  solely  on  the  joint  classification  criteria 
proposed  in  Reference  1,  in  combination  with  the  missile  Jiameters  and 
average  airframe  stiffness  properties  listed  in  the  Industry  Survey 
(Appendix  I).  The  figure  is  not  intended  to  show  actual  joint  character¬ 
istics  of  these  missiles  but  only  the  influence  that  "Excellent",  "Good", 
or  "Moderate"  joints  would  have  on  local  stiffness.  The  "Test  Specimen" 
listed  in  the  bottom  row  refers  to  the  simple  models  used  in  the  test 
program  described  in  Section  7.0. 
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Section  5,0 
PARAMETRIC  STUDIES 


5.1  JOINT  STIFFNESS  AND  LOCATION 

The  purpose  of  the  paianetric  study  discussed  in  this  section  is 
to  determine  and  illustrate  the  effects  of  variations  in  joint  stiffness 
and  location  on  missile  response  characteristics.  Emphasis  is  placed 
initially  on  studying  effects  rather  than  the  causes  of  the  various 
joint  parameters  in  order  to  arrive  at  an  understanding  of  their  rela¬ 
tive  importance. 

The  primary  tool  used  in  the  parameter  study  is  a  standard  beam 
modal  analysis  program.  This  program  accepts  a  lumped  parameter  beam 
bending  model  for  the  missile  structure  under  consideration  and  gener¬ 
ates  mode  shape,  shear,  bending  moment,  and  resonant  frequency  informa¬ 
tion.  The  approach  taken  in  determining  the  influence  of  the  joint 
parameters  on  missile  response  characteristics  has  generally  been  to 
investigate  a  single  parameter  at  a  time.  The  joint  compliance  values 
used  in  this  report  are  reciprocals  of  the  flexural  spring  rates  attri¬ 
buted  to  the  joints.  For  the  purpose  of  the  parametric  study,  only 
flexural  or  rotational  compliances  ha«e  been  considered  at  joints  since 
shear  deflections  are  generally  a  secr-.id  order  effect. 

A  non-dimensional  presentation  of  the  results  is  shown  to  provide 
a  useful  basis  for  estimating  the  effects  of  changes  in  joint  compliance 
and  location  on  the  frequencies  of  a  uniform  beam.  These  trends  devel¬ 
oped  initially  for  uniform  beams  with  single  joints  are  shown  to  be 
extrapolatible  with  surprising  accuracy  to  uniform  beams  with  multiple 
joints.  As  might  be  expected,  however,  some  discrepancies  are  likely  in 
applying  the  results  to  non-uniform  beams  representative  of  actual  tacti¬ 
cal  missiles.  One  example  based  on  ar.  actual  missile  with  five  joints 
nevertheless  indicates  that  good  qualitative  estimates  of  joint  effects 
on  the  first  mode  frequency  are  possible. 

5.1.1  Method  of  Analysis.  The  modal  analysis  program  used  in  the 
parametric  3tudy  is  based  on  a  modified  Holzer-Myklestad  method.  A 
missile  can  be  described  analytically  as  a  non-uniform  elastic  beam 
(the  main  beam)  to  which  various  appendages  are  attached.  Joint  local 
deformation  characteristics  can  be  simulated  by  placing  local  springs  at 
stations  where  joints  occur.  The  main  beam  is  represented  as  a  series 
of  lumped  masses  and  inertias  connected  by  weightless  beam  elements.  The 
beam  elements  represent  the  average  stiffness  properties  between  the  mass 
stations.  The  appendages  are  represented  in  the  same  manner  as  the  main 
beam  and  attachment  of  an  appendage  to  the  main  beam  may  be  accomplished 
with  local  springs  representing  attachment  compliances.  The  computer 
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program  is  capable  of  handling  bending,  torsion,  and  longitudinal 
motions;  however,  only  bending  analyses  for  s  free-free  beam  have  been 
considered  in  these  studies.  The  program  starts  with  an  assigned 
frequency  and  iterates  with  increasing  frequency  until  the  boundary 
conditions  are  met  and  a  natural  frequency  is  found.  At  each  natural 
frequency  found  by  the  program,  the  modal  displacement  at  some  selected 
normalization  station  (in  these  studies  the  first  station)  is  set  equal 
to  unity.  The  deflections,  slopes,  bending  roments,  and  shears  are 
obtained  at  all  missile  stations  relative  to  the  displacement  of  the 
first  station.  Hence,  all  quantities  are  obtained  as  "normalize'*" 
values . 

For  the  purposes  of  the  parametric  study  a  "nominal"  uniform  beam 
was  selected  which  has  weight,  stiffness,  and  a  first  mode  resonant 
frequency  (50.188  Hz)  similar  to  a  representative  tactical  missile. 

The  schematic  diagram  below  shows  the  "nominal"  uniform  beam  and  its 
associated  physical  properties. 


jL 


L  = 

176  IN 

M?L  = 

50.188  Kz 

D 

13.5  IS 

5.92  L3S/IN 

F.t 

2.9  (10)9  LB-IN2 

X  = 

Joint  Location  ^  IN 

This  analytical  model  approximates  the  distributed  uniform  beam  weight 
as  twenty  one  point  masses.  The  effects  of  joint  location  and  compli¬ 
ance  on  mode  shape,  bending  moment  distribution  and  on  modal  frequency 
have  been  examined  using  this  "nominal”  uniform  beam. 

5.1.2  Joint  Location  Effects.  Figure  5-1  shows  the  bending 
moment  distribution  over  the  half  span  of  the  "nominal"  uniform  beam  for 
the  first  three  bending  modes.  Only  the  half  span  distribution  is  shown 
since  the  bending  moment  is  an  even  function  about  the  midspan  station 
for  the  first  and  third  modes  and  an  odd  function  about  the  midspan 
station  for  the  second  mode.  The  bending  moment  has  been  normalized  to 
a  unity  modal  displacement  of  the  beam  at  07.  span.  Since  the  deflection 
across  a  rotational  joint  is  dependent  on  the  magnitude  of  the  bending 
moment  at  the  joint,  the  importance  of  joint  location  is  immediately 
evident.  As  shown,  the  bending  moment  reaches  a  peav  at  507.  span  in  the 
first  ar.d  third  modes  while  the  second  mode  bending  moment  gr^s  to  zero. 
Hence,  a  sort  rotational  spring  located  at  the  center  of  the  earn  wil 
produce  significant  changes  in  the  modal  characteristics  of  .ne  first 
ar.d  third  modes  while  the  second  mode  will  be  unaffected. 

Using  a  "moderate"  joint  as  determined  by  the  NASA  criteria  shown 
in  Section  4.2,  the  effect  of  joint  location  on  the  first  bending  mode 
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shape  of  the  uniform  beam  was  examined.  Figure  5-2  shows  the  first 
bending  mode  shapes  for  the  "nominal"  uniform  beam  with  a  "moderate" 
joint  at  0%  (no  joint),  10%,  20%,  30%,  40%,  and  50%  span.  Joints 
located  between  50%  and  100%  span  would  have  the  same  effects  on  the 
first  mode  characteristics  as  joints  located  between  50%  and  0%  span. 

Figure  5-3  shows  the  corresponding  first  mode  bending  moment  distri¬ 
butions  which  result  from  placing  a  "moderate"  joint  at  the  same  range  of 
locations  along  the  beam.  As  shown,  the  normalized  bending  moment 
decreases  when  a  joint  is  located  on  the  beam.  The  change  in  the  bending 
moment  distribution  is  caused  by  the  alteration  of  the  mode  shape  and 
reduction  in  modal  frequency  when  a  joint  is  added  to  the  system. 

Figure  5-4  shows  the  effect  joint  location  has  on  the  first  mode 
natural  frequency  of  the  "nominal"  uniform  beam  for  a  "moderate"  joint. 
The  maximum  reduction  in  first  mode  frequency  occurs  when  the  joint 
location  is  at  50%  span  as  might  be  expected.  A  joint  located  between 
0%  and  10%  span  causes  little  change  in  the  first  mode  frequency  while 
a  joint  location  between  40%  and  50%  span  causes  a  large  change  in  the 
first  mode  bending  frequency. 


5.1.3  Joint  Compliance  Effects.  Figure  5-5  shows  the  effect  of 
joint  compliance  on  the  first  bending  mode  shape  of  the  "nominal"  uni¬ 
form  beam  with  a  joint  located  at  50%  span.  As  expected,  the  mode  shape 
becomes  more  deformed  r.t  the  joint  location  as  the  joint  stiffness 
decreases  from  no  joint  to  a  "good"  joint  to  a  "moderate"  joint. 

Figure  5-6  shows  the  resulting  bending  moment  distributions  for  no  joint, 
for  a  "good"  joint,  and  for  a  ’•moderate"  joint.  Again  the  distribution 
becomes  increasingly  altered  as  the  joint  spring  becomes  softer. 

Figure  5-7  shows  the  effect  of  joint  stiffness  on  the  reduction  of  the 
first  mode  bending  frequency  due  to  joint  location. 


5.1.4  Non -dimensional  Presentation.  Based  on  the  uniform  br  n 
study,  a  useful  non-dimensional  relationship  can  be  derived  betweu.i 
joint  compliance  and  beam  stiffness.  This  relationship  provides  a  means 
of  easily  assessing  the  significance  of  altering  a  joint  parameter.  The 
slope  change  across  a  rotational  joint  is  equal  to  the  product  of  the 
bending  moment  at  the  joint  and  the  joint  compliance  {  4^/  MC,^  ). 

The  change  in  slope  per  unit  length  on  a  uniform  beam  without  joints  is 
equal  to  the  bending  moment  divided  by  the  beam  stiffness 
A  compliance  ratio  per  unit  beam  length  is  obtained  by  dividing  the  joint 
rotation  per  unit  moment  by  the  beam  rotation  per  unit  length, 
a<4,'A«k;  m  *=  .  This  term  can  be  conveniently  non -dimens  ional- 

ized  by  dividing  by  the  total  beam  length  JL  .  A  nondimensional  compliance 
ratio, >  is  thus  obtained  which  relates  joint  compliance,  beam 
stiffness £  and  bean  length.  This  expression,  not  surprisingly,  is 
directly  related  to  the  stiffness  reduction  ratio  derived  in  Section  4.1. 
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Using  total  beam  length  to  non -dimens ionalize  the  compliance  ratio 
offers  the  advantage  of  general  application  to  all  uniform  beams  regard¬ 
less  of  fineness  ratio  (length  divided  by  diameter).  The  use  of  body 
diameter  as  a  basis  for  non-dimensionalizing,  however,  is  believed  to 
offer  better  physical  insight  in  judging  the  significance  of  joint  com¬ 
pliance  relative  to  airframe  stiffness.  A  reference  length  baaed  on 
diameter  also  permits  direct  comparison  with  the  joint  compliance  rating 
system  discussed  in  Section  4.2.  The  results  presented  on  this  basis, 
however,  apply  explicitly  only  for  airframes  with  a  fineness  ratio  of 
13  as  used  in  the  parametric  study.  Both  non-dimenslonalizing  parameters 
are  given  in  Tables  5-1,  5-2  and  5-3  which  summarize  in  tabular  form  the 
effects  of  non-dimensional  joint  compliance  and  location  on  the  mode  fre¬ 
quencies  of  a  uniform  beam.  Corresponding  plots  of  the  first  mode 
frequency  ratio  versus  compliance  ratio  (EiC/L)  and  stiffness  ratio 
(Kr)  are  presented  in  Figures  5-8  and  5-9  respectively.  It  is  signifi¬ 
cant  to  note  that  a  single  "moderate*'  joint  at  the  midspan  of  a  uniform 
beam  can  be  expected  to  reduce  the  first  mode  frequency  by  approximately 
35  percent. 

5.1.5  Multiple  Joint  Applications.  The  uniform  beam  results  for  a 
single  joint  as  displayed  in  Figures  5-8  and  5-9  can  be  extrapolated  to 
multiple  joint  configurations  by  a  simple  summing  procedure  in  which  a 
number  of  joints  of  various  compliances  at  arbitrary  locations  are  refer¬ 
enced  and  summed  as  an  equivalent  single  joint  compliance  located  at  a 
common  reference  station  such  as  the  midspan.  The  steps  in  this  proced¬ 
ure  are  as  follows 

1.  At  each  desired  joint  location  determine  the  resulting  frequency 
ratic  for  a  single  joint  of  specified  compliance  at  that 
location. 

..  Determine  the  equivalent  joint  compliance  which  would  be 

required  at  a  common  reference  station  such  as  the  midspan  to 
produce  the  same  frequency  ratio. 

3.  The  sum  of  the  individual  equivalent  joint  compliances  yields 
an  equivalent  single  joint  at  the  reference  station. 

4.  The  frequency  ratio  associated  with  this  single  equivalent 
joint  is  found  to  give  a  very  close  estimate  of  the  multiple 
joint  effect  on  the  first  mode  frequency  ratio  for  a  uniform 
beam. 


An  example  of  the  application  of  this  procedure  is  presented  in 
Table  5-4.  The  "calculated"  frequency  ratios  for  the  multiple  joint 
configurations  listed  were  determined  by  modal  analysis,  with  the 
"estimated"  values  determined  by  the  procedure  described.  In  this 
instance,  the  use  of  compliance  ratio,  EIC/L,  offers  an  advantage  over 
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stiffness  ratio,  KR,  in  that  the  incremental  equivalent  compliances  can 
be  added  directly. 

5.1.6  Tactical  Missile  Application.  For  an  actual  tactical  missile, 
non-uniformities  exist  in  both  the  mast?  and  stiffness  distribution  in 
addition  to  multiple  joints.  The  effect  of  these  considerations  would  be 
expected  to  severely  limit  the  applicability  of  the  uniform  beam  results 
in  estimating  joint  effects  on  an  actual  missile.  The  following  example, 
however,  illustrates  that  a  good  qualitative  estimate  of  joint  compli¬ 
ance  effects  is  nevertheless  possible.  Figure  5-10  displays  a  typical 
a^r^rame  stiffness  distribution,  in  this  case  for  the  medium  range 
Standard  Missile.  The  significance  of  the  narrow  El  spikes  showing  dis¬ 
continuous  increases  in  airframe  stiffness  which  are  often  meticulous lv 
computed  and  included  in  modal  analysis  computations  is  virtually  zero, 
having  very  little  effect  in  theory  as  shown  by  the  calculated  frequency 
changes  in  the  figure.  This  conclusion,  however,  definitely  does  not 
apply  in  the  case  of  discontinuous  decreases  in  El  as  will  be  demonstrated 
in  the  test  results  presented  in  Section  7.  The  following  table  lists 
the  five  principal  airframe  joints,  their  compliances  and  locations,  and 
the  comparison  between  the  computed  (modal  analysis)  and  estimated  fre¬ 
quency  ratios  by  the  procedure  of  Section  5.1.5. 


STANDARD  MISFILE  (MR)  JOINT  CHARACTERISTICS 


Joint 

Location 

Compliance 

Ce  x  (10)8 

Compliance  Ratio  EIC/L(10)2 

Actual 

Equiv. 

.22 

1.5 

24.70 

5.06 

.35 

0.6 

9.85 

6.52 

.41 

0.5 

8.25 

6.90 

.87 

1.15 

18.95 

0.88 

.92 

1.0 

16.45 

0.31 

j  1st  Mode  Frequency  Ratio 

Error 

7. 

Calculated 

Estimated 

.834 

.810 

3 

ihe  implied  accuracy  of  the  estimate  is  somewhat  misleading  since 
an  average  airframe  stiffness  must  be  assumed  for  the  total  airframe  and 
the  answer  is  obviously  .ensitive  to  this  assumption.  The  answer  does 
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suggest,  however,  that  the  procedure  may  be  useful  in  providing  good 
qualitative  estimates  of  joint  effects. 

5.2  CROSS  PLANE  SPRING  COUPLING  AT  AIRFRAME  JOINTS 

5.2.1  Introduction.  Double  peaked  modes  are  cften  observed  in  the 
vibration  testing  of  tactical  missiles.  This  condition  is  most  usually 
accompanied  by  substantial  energy  transfer  or  cross  talk  between  in- 
plane  and  cut-of-olane  response.  The  cause  of  this  coupling  can  frequent¬ 
ly  be  traced  to  the  behavior  of  one  or  more  of  the  airframe  joints. 
"Breaking*1  and  retorqueing  of  the  joints  will  in  many  instances  signifi¬ 
cantly  alter  or  in  some  cases  eliminate  this  response  characteristic. 

One  of  the  major  problems  created  by  this  type  of  mechanical  joint  behav¬ 
ior  is  the  uncertainty  it  introduces  in  describing  missile  airframe 
response  for  control  system  design  studies.  One  typical  frequency 
response,  measured  on  the  nose  of  a  grain-out  extended  range  Standard 
Missile  for  constant  force  excitation  at  the  tail  station,  is  prest  ;ed 
in  Figure  5-11.  Peaks  occur  at  both  63  Hz  and  66  Hz.  The  peak  at  the 
lower  frequency  is  1.5  db  or  about  20  percent  greater  than  the  peak  at 
the  higher  frequency.  It  has  been  determined  that  for  this  airframe  the 
severity  of  the  double  peak  characteristic  is  primarily  controlled  by  the 
forward-most  joint  w’'ich  is  a  split  ring  continuous  land  design  (Figure 
1-15).  Reseating  and  overtorqueing  of  the  coupling  nut  will  generally 
bring  the  peaks  closer  together  in  frequency  and  increase  the  magnitude 
of  the  peak  response. 

The  objective  in  this  section  is  to  develop  a  conceptual  model  for 
the  source  of  elastic  coupling  in  airframe  joint3  and  to  examine  joint 
coupling  effects  on  airframe  response  through  parameteric  analysis. 

5.2.2  Conceptual  Model  for  Joint  Elastic  Coupling.  It  is  believed 
that  joint  elastic  coupling  is  caused  by  irregularities  in  the  mating 
surfaces  of  the  joint.  These  irregularities  produce  a  non -axially 
symmetric  load  distribution  in  the  joint  and,  hence,  create  spring  coup¬ 
ling  characteristics.  Discrete  irregularities  can  be  idealized  as 
springs  concentrated  at  the  points  of  protrusions  around  the  airframe 
circumference.  If  concentrated  springs  (with  spring  rates  k^)  are 
assumed  to  be  located  in  the  plane  of  an  airframe  joint  as  shown  below, 
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an  equation  can  be  derived  expressing  the  bending  moment  across  the 
joint  in  the  ’Z1  and  'Y*  directions  as  a  function  of  the  spring  rates 
,  angles  ,  and  joint  rotations  ot ^  and  ocy  .  The  expression  for 
bending  moment  across  a  joint  is 


Z  K-  SM/* 

-Z  *£  3M/&S  <3&sa>t. 


-Z  a;-  s/A/&i  cos&. 
Z  K  c.esx&i 


M  =  bending  moment  across  the  Joint 

K  =  rotational  stiffness  of  the  joint 

=  rotation  across  the  joint 

3E,Y-  two  reference  missile  axes  which  are  perpen¬ 
dicular. 


A  useful  relationship  for  investigating  joint  coupling  effects  can 
be  developed  by  reducing  the  general  case  to  a  three  spring  idealization. 
Since  three  points  define  a  plane,  this  idealization  offers  a  plausible 
model  for  the  points  of  joint  interface  contact  and  becomes  especially 
valid  when  the  source  of  the  preload  in  the  joint  (e.g.  bolts  in  tension) 
represents  a  minor  contribution  to  joint  compliance. 


Considering  the  plane  of  the  joint  for  the  three  spring  idealization, 
and  equal  spring  rates: 


K  -  *3 


POMONA  OMRATtON 


OMNBRAL  DYNAMIC* 

Etactro  Dynamic  Division 


The  spring  matrix  in  the  y*.  coordinate  system  becomes: 


where:  , 

*  =  T  (*'***  +  '* 


ana  ^  v 

aK  -  C  /  ~  “  SMS  «S>)  K 

It  is  significant  to  note  that  for  all  values  of  0  a  plane  of  symme¬ 
try  exists  (defined  by  the  Y  axis  and  the  missile  longitudinal  axis) 
and  that  the  Y-X.  coordinates  remain  decoupled.  The  special  case  of 
axial  symmetry  (  aK*0)  exists  when  0  =  60°  (120°  spacing  between 
springs) . 

The  joint  deflection  characteristics  viewed  from  a  y'~m‘  coordinate 
reference  rotated  through  45  degrees  will  exhibit  maximum  apparent  coup¬ 
ling  through  -»/C  *n  the  following  transformed  spring  matrix: 


Figure  5-12  presents  the  variation  in  with  spacing  angle  0 

for  the  three  spring  idealization  joint  model.  It  should  be  noted  thct 
the  limiting  case  of  100  percent  coupling,  aX*-ac  implies  a  zero  spring 
about  the  oty  coordinate  and  maximum  coupling  in  the  y-  m *  coordinate 
frame. 

For  actual  airframes,  particularly  with  multiple  joints,  coupling 
effects  on  airframe  dynamic  response  car.  become  sizable  and  quite  com¬ 
plex.  Although  each  joint  may  have  a  plane  of  symmetry  with  respect  to 
its  load /defe let ion  characteristics,  the  probability  of  alignment  in 
multiple  joint  planes  of  symmetry  is  quite  small. 
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5.2.3  Joint  Spring  Coupling  Parametric  Study.  A  computer  program 
which  accepts  generalized  modal  properties  for  a  linear  damped  system 
and  computes  the  frequency  response  of  the  desired  structural  transfer 
function  has  been  used  in  the  parametric  study  of  joint  spring  coupling. 
The  spring  coupling  term  governs  the  response  in  any  one  plane  which  is 
solely  due  to  an  input  in  a  plane  perpendicular  to  the  plane  of  response. 

The-  equations  of  motion  of  a  linear  system  undergoing  forced 
vibrations  may  be  written  in  matrix  form  and  Laplace  operator  notation 
as  follows: 


where 


K]  **  *[*]«”■  M]  ft-3  *  frJ 


[*] 

[«>■] 

M 


generalized  mass  matrix 


=  generalized  damping  matrix 


=  generalized  spring  matrix 


=  time  varying  generalized  coordinates  of  the  system 


generalized  forces  acting  on  the  system 


dx 

the  Laplace  operator,  SX  •  . 


If  ^  is  the  time  varying  generalized  coordinate  representing  the 
forcing  function. 


tel  -  H  K, 


and  substituting. 


|t<] s* *  m  j  *  w  *s- 


The  response  of  the  system  is  a  linear  combination  of  the  coordinate 
responses,  which  is  found  by  Cramer’s  Rule: 
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Sf 

where 

0(s> 

- 

the  determinant  expansion  of  the  left-hand  side  of 
the  equations  of  motion 

a 

the  determinant  expansion  of  the  left-hand  side  of 
the  equations  of  motion  after  substituting  the  right- 
hand  side  into  the  1th  column 

= 

a  numerical  coefficient  (modal  deflection  at  point 
of  interest  on  beam) 

X 

= 

the  sectional  response  of  the  system  to  the  forcing 
function  3  f.  . 

If  the  response  is  desired  in  terms  of  velocity  or  acceleration, 
the  response  equation  may  be  multiplied  by  s  or  s2,  respectively: 


9f  =  DCs)  »  sf  *  *>C*> 

The  generalized  coordinates  used  in  the  parametric  study  consist  of 
in-plane  and  out-of -plane  bending  modes  and  in-plane  and  out -of -plane 
"hinge''  modes.  The  bending  modes  have  no  compliance  at  the  Joint  of 
Interest,  and  consequently  have  no  change  of  slope  across  the  joint. 

The  "hinge"  modes  treat  the  beam  as  two  rigid  members  connected  with  a 
local  spring  at  the  joint  of  Interest.  The  local  spring  represents  the 
joint  compliance  in  the  two  individual  planes  and  the  joint  coupling 
characteristics  between  the  planes. 

The  displacaments  in  the  "hinge"  modes  are  such  that  they  are  mass 
orthogonal  with  respect  to  rigid  body  translation  and  pitch  about  the 
center  of  gravity.  The  "hinge"  modes  are  not  orthogonal  to  the  bending 
modes,  however,  and  hence  in  each  plane  we  have  mass  coupling  terms 

*  MJ£  *  f m 

o 


where 

m  =  running  mass  (mass  per  unit  length) 
*  ith  bending  mode  shape 
-  "hinge"  mode  shape 
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The  only  coupling  between  the  two  planes  of  motion  occurs  solely  in 
the  structural  spring  matrix.  When  Lagranges1  equations  are  used  to 
derive  the  equations  of  motion,  it  can  be  seen  that  the  structural 
spring  terms  result  from  the  stored  energy  associated  with  the  unit  deflec¬ 
tions  of  all  the  generalized  coordinates.  The  hinge  modes  store  energy 
only  by  rotation  at  the  joint,  and  since  the  bending  modes  have  no  change 
of  slope  at  the  joint,  there  is  no  spring  coupling  between  the  bending 
and  hinge  modes.  Consequently,  Che  spring  matrix  for  the  bending  modes 
is  a  diagonal  matrix.  The  two  hinge  modes  have  both  diagonal  terms  and 
off  diagonal  terms,  which  represent  energy  transfer  between  planes.  The 
spring  terms  for  the  hinge  modes  are  given  by  the  following  relations: 


&U  C  *aa  =  (+  <t>-y  Y  Xyy 


where 


=  change  of  slope  across  the  joint  for  the  hinge 
modes  A  *Jy  ) 


Kzz  =  j°int  spring  in  the  Z  plane  (Kyy  =  Kzz) 

Kzy  =  chosen  as  a  parameter  (a  percentage  of  Kzz) 

i,j  =  mode  counters  of  hinge  mode  in  z  and  y  planes, 
respectively. 


The  damping  matrix  used  in  this  study,  M  >  !■  a  diagonal  matrix. 
The  and  QqJ  matrices  have  the  following  form  when  considering 

three  bending  modes  and  one  "hinge"  mode  in  both  the  y  and  z  planes: 


X 

X 

X 

XXX 


X 


X  X 


X 

X 
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I 
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J 


where 


First  bending  mode  in  Z -plane 
Second  bending  mode  in  Z-plane 
Third  bending  mode  in  Z-plane 
"Hinge”  mode  in  Z-plane 
First  beading  mode  in  Y -plane 
Second  bending  mode  in  Y-plane 


Third  bending  mode  in  Y-plane 
=  "Hinge”  mode  in  Y-plane 


very  convenient ^or^nar^V*1?0  Set  of  generalized  coordinates  is 

Any  change  of  joint  coJoU™  n  ?  ““Vf  J°lnt  ccu',11"8  ^oractcristicn. 
than  couple  thi  hLge  “=Pa  5”S  uht  rV^  tems  “f  the  KI  »‘rlx 
compliance  in  the  planes  or’tn^^  iL<eSired  chan8e  the  Joint 
this  again  only  involves  the  spring  tenS^VH* dances  in  the  planes, 
Consequently,  the  effect  of  pertatninS  t0  the  hinge  modes, 

investigated  with  fminLm  If  f  °f  the  j°int  Pities  can  be 
of  motion.  nlmUra  °f  Change  to  the  overall  set  of  equations 


acteristicSfofta  wif cmTbeam^8  C°?pling  on  the  dyna.i£c  response  char- 

"nominal"  beam  of  sJctl™  5  T  3"?.,  ? tu?led  Analytically  „si„8  the 

section  5.1.3  with  a  single  local  flexural  compliance. 
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The  bending  frequencies  and  mode  shapes  for  these  studies  were  computed 
using  the  modified  Holzer-Myklestad  computer  program  discussed  in 
Section  5.1.1.  A  damping  ratio  of  one  percent  of  critical  was  assumed 
for  all  modes  and  the  frequency  response  curves  were  computed  and  plotted 
using  the  transfer  function  program  described  earlier  in  this  section. 

The  amount  of  cros.  plane  spring  coupling  is  denoted  as  a  percentage  of 
the  in-plane  joint  spring  term. 

Figure  5-13  shows  the  effect  of  spring  coupling  on  the  frequency 
response  characteristics  of  a  uniform  beam  with  a  ''good"  joint  at  50% 
span.  The  acceleration  response  was  computed  for  a  point  at  zero  span 
for  force  excitation  at  100%  span.  The  various  curves  are  for  zero, 

50%  and  75%  cross  plane  spring  coupling.  As  can  be  seen,  the  peak 
response  decreases  with  increasing  spring  coupling  and  doubled  peaked 
modes  occur.  Increasing  attenuation  exists  at  the  center  frequency  as 
the  spring  coupling  increases  (the  valley  between  the  two  peaks  becomes 
deeper).  Figure  5-14  shows  the  effect  of  joint  spring  coupling  on  the 
response  characteristics  of  a  uniform  beam  with  a  "moderate"  joint  at 
50%  span.  The  effects  of  spring  coupling  are  much  more  pronounced  for 
a  "moderate"  joint  than  for  a  "good"  joii t.  Since  tie  bending  moment  in 
the  second  mode  is  zero  at  the  joint  location  for  both  of  these  cases  no 
effects  of  joint  stiffness  or  coupling  are  seen  in  the  second  bending 
mode . 

The  trend  in  the  magnitude  of  joint  coupling  effect  on  the  first 
bending  mode  frequency  as  the  joint  compliance  is  varied  is  illustrated 
in  Figure  5-15.  In  this  figure,  the  total  frequency  shift  ) 

as  a  function  of  joint  elastic  coupling  is  shown  for  both  the  "moderate" 
and  '  goed"  joint  located  at  the  midspan  of  the  nominal  uniform  beam. 

The  point  to  be  made  is  that  relatively  small  amounts  of  elastic 
coupling  in  mec  lanical  joints  can  produce  substantial  changes  in  both 
gain  and  resonant  frequency  characteristics.  Conti’ol  system  body  mode 
coupling  filters  typicall /  designeu  for  a  body  mode  frequency  tolerance 
of  ±5  percent  may  be  forced  to  cope  with  a  much  larger  range  of  response 
uncertainties . 

5.2.4  Spring  Coupling  Studies  ior  an  Actual  Tactical  Missile, 

As  was  pointed  out  earlier,  Figure  5-11  shows  a  measured  Standard  Missile 
ER  frequency  response  for  a  force  excitation  at  the  tail  station,  A 
study  has  been  conducted  which  attempts  the  modeling  of  the  missile  in 
such  a  manner  as  to  predict  a  response  similar  to  that  of  Figure  5-11. 

The  stiffness  coupling  terms  were  calculated  usfng  two  sets  of  assumed 
modes  which  represent  the  missile  dynamic  characteristics  in  two  ortho¬ 
gonal  planes.  The  two  sets  consisted  of  four  modes  each.  The  frequency 
response  of  the  resulting  equation*;  of  motion  was  computed  using  the 
computer  program  previously  mentioned.  The  2-plane  frequency  responses 
to  2 -plane  force  input,  calculated  /ith  and  without  spring  coupling 
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|  between  the  ,JZ"  and  "Y"  planes,  are  shown  in  Figure  5-16.  The  case  in 

|  which  no  coupling  exists  between  the  two  planes  shows  a  single  sharp 

|  peak  in  the  first  mode.  The  case  in  which  elastic  coupling  is  introduced 

|  between  the  Z-plane  and  the  Y -plane  equal  to  50  oercent  of  the  total 

|  generalized  spring  is  shown  to  provide  a  double  peaked  response  in  close 

1  agreement  with  the  measured  data  of  Figure  5-11.  It  is  felt  that  this 

1  joint  stiffness  coupling  is  a  first  order  model  of  an  important  deviation 

I  from  the  single  plane  response  model  used  in  classical  structural  dynamic 

1  analysis  of  tactical  missiles . 
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Table  5-1 

FIRST  MODE  FREQUi  .'ICY  RATIO  VARIATION 
WITH  JOINT  LOCATION’  AND  RATING 


8  B 

■s 

Joint  Location 

mraQH|| 

0.2 

0.3 

BTrmr 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

£ 

.883 

1.0 

.999 

.997 

.994 

,993 

.  ,599 

1.0 

.997 

.990 

.983 

.979 

G 

.429 

.999 

.990 

,970 

.949 

.940 

.189 

.997 

.969 

.909 

.857 

.838 

IV 

.0698 

.991 

.904 

.768 

.682 

.655 

.0227 

.970 

.730 

.541 

.457 

.433 

L 

:  .00745 

.891 

.490 

.337 

,279 

.263 

.00232  t 

.648 

*  " 

.289 

.194 

.159 

.150 

Table  5-2 

MODAL  FREQUENCY  RATIO  VARIATION 
WITH  MIDSPAN  JOINT  RATING 


Table  5-3  - 

MODAL.  FREQUENCY  RATIO  VARIATION 
WITH  "MODERATE"  JOINT  -LOG \TION 


T-*bl«  5-4 

UNIFORM  BEAM  WITH  MULTIPLE  JOINTS 
CALCULATED  AND  ESTIMATED  FREQUENCY  RATIOS 
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Section  6.0 

EXTRACTION  0?  JOINT  COMPLIANCES  FROM  MEASURED  MODAL  DATA 


Computed  modes  and  frequencies  of  a  tactical  missile  containing 
joints  often  will  vary  considerably  from  experimental  mode  shapes  and 
frequencies.  Since  the  mass  and  basic  stiffness  properties  of  the 
missile  are  known  with  a  good  degree  of  certainty,  the  values  used  for 
the  joint  compliances  must  account  for  the  discrepancies  in  most  cases. 
If  there  are  multiple  joints  and  several  modes  involved,  the  task  of 
iteratively  adjusting  compliance  values  by  hand  to  produce  matching 
results  becomes  laborious,  time  consuming,  and  often  very  frustrating. 
Consequently,  It  is  desirable  to  automate  a  method  that  will  extract 
the  joint  compliances  in  a  manner  such  that  theoretical  results  match 
experimental  results  as  closely  as  possible.  Ihe  following  section 
of  this  report  presents  the  equations,  evaluation,  and  some  limitations 
of  an  initial  effort  to  develop  such  a  capability. 

6.1  METHOD  OF  ANALYSIS 

A  method  currently  being  used  to  compute  optimum  trajectories  is 
applicable  to  the  problem  (References  2  and  3).  It  is  called  the  method 
of  steepest  ascent,  and  it  uses  the  steepest  linear  slope  in  a  function 
space  to  proceed  towards  the  minimum  (or  maximum)  of  seme  computed 
quantity.  As  applied  to  the  problem  of  adjusting  joint  compliances, 
it  searches  for  a  distribution  of  control  variables  (joint  compliances) 
and  unspecified  initial  conditions  (modal  slopes)  that  will  give  an 
optimum  for  a  quantity  (quadratic  function  of  mode  shape  errors)  that  is 
a  function  of  the  integrated  stats  variables  (modal  deflections)  sub¬ 
ject  to  constraints  (boundary  conditions). 

Ih  general  form,  the  matrix  differential  equations  to  be  solved 
for  beam  modal  deflections,  internal  forces,  constraint  quantities,  and 
optimizing  quantity,  can  be  expressed  as 


fw/}  « (g  C  k-  .  »  *>}  *  <Vi,<*0*)]  £*}  4  {k  C«<  .*>} 


^  s  * 

where 

=  computed  (state)  quantities 

c  y 

=  differentiation  with  respect  to  X 

3  control  variables  (quantities  to  be  adjusted) 

M, 

ft} 

-  coefficient  matrices 

X 

=  coordinate  along  the  beam 
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If  estimates  are  made  for  the  initial  values  of  the  state  variables  (w/) 
and  for  values  of  the  distributed  control  variables  (at),  the  equations 
can  be  integrated  numerically  along  the  length  of  the  missile.  In 
general  the  boundary  conditions  will  not  be  met  and  the  error  ^optimize- 
tion)  function  will  not  be  a  minimum.  It  is  required  to  alter  the  con¬ 
trol  variables  («),  in  a  method  chat  tends  to  satisfy  the  desired  end 
results.. 


For  small  change?  **  ■  the  Hoft-tr  kltC  previous 

solution  can  be  described  by  the  linear'  perturbation  equations 

■-  fa  fa  ,«<>*)]  (fa}  *  | %  fa  >  .  *)]  (s«i} 

where 


The  set  of  equations  adjoint  to  the  perturbation  equations  are  defined 
as 


{*•'}  -  -  EWfr 

M  -  adjoint  solutions 
[  ]  3  matrix  transpose 


These  adjoint  solutions  are  related  to  the  perturbation  solutions  by 
the  relationship 

(L\J  f<CL*J[f<j3£r«i})  *  *  CL*J 

where 


«  =  OtJi  3  beginning  and  end  of  the  missile. 


This  relationship  proves  very  useful  for  properly  selected  solutions  of 
the  adjoint  equations.  Lee  the  quantity  of  interest  at  J2  be  designated 

If  we  define  the  value  of  A,-  Ci)  as  — ,  the 

3  faOO) 
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adjoint  equations  can  be  integrated  backward  from  XmJf-  to  jt  to 

give  solutions  for  M  •  The  left  hand  side  of  the  equation  relating 
the  perturbations  to  the  adjoint  solution  can  be  seen  to  equal  the 
perturbation  of  for  this  set  of  JV’a  .  Looking  at  the  right  hand  side 
of  the  equation,  the  are  influence  coefficients  that  relate  changes 

in  the  state  variables  at  x  »-<S>  and  control  variables  to  changes  in  $  . 
The  method  of  steepest  ascent  uses  this  relationship  to  choose  the  manner 
to  vary  <**•  and  Co )  in  order  to  proceed  rapidly  to  the  desired  value 

of  #  . 


The  particular  steepest  ascent  method  followed  is  the  second  method 
outlined  in  Reference  3;  Dropping  matrix  notation,  a  generalized  weighted 
mean  square  "step  size"  is  defined  as 

/ 

C V/O*  ~  Am*  M* 2  :  +  J 0*04 <*)  dx 

c 

where 


dP  -  step  size 

Am  =  a  column  of  changes  to  hj/s) 

(only  involves  those  that  are  variable) 

rf'oc  =  running  change  in  the  control  variables 

j£tO  *-  weighting  function  matrices  - 

C  /  *  matrix  transpose 

The  desired  amount  of  change,  to  the  payoff  quantity  (+$  )  and  the 
necessary  changes  to  meet  the  constraints are  specified  and  the 
(4P)*  necessary  to  accomplish  this  is  minimized.  This  results  in  the 
following 


*  * 


i'*t'  r~!. 

*/■»■  me*- 


o'  o'  r*~'  a-t_ 

*/*■  */a 


where 


w~'  gt~  ' ^ 
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V-  ‘  [\>  I  *  J 


(influence  functions  corresponding 
to  changes  in  #  and  to  the 
constraints ,  -*•*  ) 


-  [*4»  ! 


(only  those  corresponding  to  w^Co) 
that  are  variable) 


+  *»  (previously  defined) 

The  a*  are  used  to  alter  the  initial  conditions  on  the  state 
variables ,  C^i  ^  »  and  used  to  alter  the  control  variables  along 

the  missile.  The  differential  equations  are  then  integrated  again  to 
give  a  new  solution  for  and  the  process  repeated  until  a  minimum 

is  reached  for  £  ,  the  payoff  quantity.  If  the  differential  equations, 
$  ,  and  JL  are  all  linear  with  respect  to  the  state  variables  it  only 
takes  one  pass  to  reach  the  answer.  (In  such  a  case,  f  is  actually  a 
constraint,  not  an  optimizing  function). 

This  general  method  was  applied  to  the  problem  of  finding  joint 
compliance  values  from  beam  bending  mode  data.  The  state  variables 
have  the  form 


where  each  submatrix  corresponds  to  an  input  mode.  For  each  mode  or 
submatrix  we  have 


M"  - 


W~~  Deflection 


W ~  Slope 
BM~  Bending  Moment 


S  Shear 


nth  mode  payoff  quantity 
7*—  Mass  coupling  to  translation 
Mass  coupling  to  pitch 
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Tl.e  only  type  of  control  variable  that  we  have  is  joint  compliance, 
consequently 


=  Cg,  a  scalar 

where  Cg  =  Bending  compliance 

Two  different  sets  of  differential  equations  were  used.  The  first 
uses  input  mode  shapes  to  describe  the  loadings  on  the  missile,  and  the 
second  uses  the  computed  shapes.  In  both  cases  a  lumped  mass  model  with 
constant  stiffness  between  masses  was  used.  The  lumped  masses  and  the 
concentrated  joint  compliances  result  in  discontinuities  in  some  of  the 
differential  equations,  but  this  does  not  cause  a  problem  with  numerical 
integration. 

For  the  first  set  of  equations,  we  have  (for  each  mode) 
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where 

=*  bending  stiffness 
<t>n  *  nth  mode  shape 

C  ^  “  quantity  at  a  station  k 

*n  *  mass 

X  *  distance  from  center  of  gravity 

&?n  m  frequency  of  nth  mode 

The  second  set  of  equations  varies  from  the  first  only  in  that 

(**,Y  -  -  »»«  <*>n 

and 

<?>)"  -  o- 


The  boundary  values 


W'-  for  the  Initial  pass  are 
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The  only  boundary  values  that  are  allowed  to  vary  are  the  a9  which 
are  adjusted  by  ir.  subsequent  iterations.  The  total  payoff  quantity 
is  taken  to  be 


#  -  r  a„\ 

and  the  constraints  (one  for  each  mode)  are 


*  A 


The  adjoint  equations  are  integrated  backwards  to  obtain  the  influence 
function,  once  for  those  related  to  the  payoff  quantity  and  once  for 
each  of  the  constraints.  The  initial  conditions  for  those  integratlcns 
are  defined  by 


for  ths  payoff  function  and 


for  each  of  the  constraints.  The  influence  functions  are  then  used  to 
calculate  changes  to  the  initial  deflection  slopes  and  to  the  bending 
compliances  in  a  manner  that  will  reduce  $  and  tend  to  match  the 
constraints.  The  process  is  repeated  until  #  is  minimized  to  within  a 
specified  tolerance. 
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Both  formulations  of  the  deflection  equations  were  used  in  demon¬ 
strating  the  feasibility  of  the  overall  method.  The  first  formulation 
has  the  advantage  that  the  input  mode  shapes  are  reflected  in  the  beam 
loadings  as  well  as  the  payoff  quantity.  It  was  used  to  make  initial 
alterations  to  the  compliances.  In  the  second  formulation  the  input  mode 
shapes  appear  only  in  the  payoff  quantity.  The  deflections  are  computed 
in  the  same  manner  as  would  be  done  in  a  beam  modal  analysis  and  conse¬ 
quently  represent  computed  modes.  In  both  cases  the  measured  frequen¬ 
cies  are  used  and  are  considered  not  to  be  variable.  The  frequencies 
could  be  treated  as  control  variables  with  a  relatively  simple  altera  - 
tion. 

6.2  TEST  CASE 

A  test  case  was  developed  for  the  purpose  of  exercising  and  checking 
out  the  joint  compliance  extraction  technique.  An  analytical  “test" 
case  was  developed  which  was-  free  of  experimental  error.  The  figure  and 
table  below  show  this  test  case  which  consists  of  an  idealized  missile 
with  four  uniform  and  one  nonuniform  beam  sections  connected  by  four 
local  flexural  springs  representing  four  joints. 
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Figure  6-1  3bows  the  weight  and  stiffness  distributions  selected  for  the 
test  case.  A  lumped  parameter  representation  as  discussed  in  Section 
5.1  was  used  to  compute  the  first  four  bending  mode  shapes  and  natural 
frequencies  for  use  in  the  development  of  the  compliance  extraction 
technique.  The  test  model  natural  frequencies  with  and  without  the  four 
Joints  are: 
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Natural  Frequencies 
With  Joints 
(Hz) 

Natural  Frequencies 
Witout  Joints 
(Hz) 

Frequency 

Ratio 

40.3 

50.4 

.800 

104.5 

118.3 

.883 

196.1 

231.5 

.847 

313.3 

389.9 

.804 

It  can  be  seen,  from  the  frequency  changes  shown  above,  that  the  joints 
selected  for  the  test  case  have  a  significant  influence  on  the  beam  bend 
Lng  modes  and  are  reasonably  typical  of  actual  tactical  missiles. 

The  first  four  bending  mode  shapes  and  frequencies  were  utilized  in 
the  joint  compliance  extraction  computer  .program.  The  stiffness  and 
weight  distributions  were  assumed  to  be  known  correctly.  Erroneous 
initial  values  were  assumed  for  the  four  joint  compliances,  with  the 
maximum  error  being  a  factor  of  ten  from  the  actual  compliance  values. 
The  only  unknowns  were  the  four  joint  compliance  values.  Figure  6-2 
shows  the  error  in  the  originally  assumed  compliance  of  the  four  joints, 
'  id  the  error  which  existed  after  each  of  the  six  computational  itera¬ 
tions  .  The  first  three  iterations  used  the  input  mode  shapes  in  solving 
for  the  four  joint  compliances.  The  second  three  iterations  used  the 
mode  shapes  generated  by  the  compliance  extraction  computer  program  in 
solving  for  the  compliance  values.  The  first  procedure  produces  better 
results  than  the  second  procedure  for  this  .particular  test  case  because 
a  perfect  set  of  input  mode  shapes  are  being  used.  In  general  this  will 
not  be  the  case  due  to  the  experimental  error  in  the  mode  shapes.  As 
shown,  the  convergence  of  the  solution  for  the  joint  compliance  values 
Is  very  rapid.-  Only  cue  Iteration  is  required  for  the -program  to  con¬ 
verge  on  the  solution.  The  computed  joint  compliance  with  the  largest 
error  is  that  of  the  second  joint  which  also  has  the  smallest  compliance 
of  the  four  joints.. 

6.3  SENSITIVITY  ANALYSIS 

The  joint  compliance  extraction  program  is  seen  to  be  functioning 
properly  for  the  test  case  discussed  above.  However,  real  test  data 
will  vary  from  this  ideal  test  case  because  of  the  presence  of  experlmen 
tai  error  in  the  data.  Errors  will  be  present  in  both  the  measured  mode 
shapes  and  frequencies.  The  sensitivities  for  the  joint  compliance  of 
the  test  model  to  certain  types  of  errors  In  the  input  data  have  been 
investigated. 

Table  6-1  shows  the  effects  of  input  modal  frequency  variations  on 
the  four  joint  compliances  of  che  test  model.  The  frequencies  of  all 
four  modes  were  varied  simultaneously  by  amounts  of  plus  1%,  plus  21, 
minus  11,  and  minus  27..  The  four  joint  compliances  have  a  moderate 
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dependence  on  the  input  frequencies. 

The  sensitivities  of  the  four  joint  compliances  to  the  number  of 
input  modes  are  shown  in  Table  6-2.  At  the  present  time,  the  results 
are  highly  dependent  on  the  number  of  measured  modes.  The  zero  values 
in  Table  6-2  are  the  results  of  the  program  trying  to  assume  a  negative 
value  for  these  particular  compliances.  Four  joint  compliances  are 
present  in  the  te3t  case,  and  four  input  modes  are  required  to  produce 
good  results.  If  the  number  of  modes  considered  is  less  than  the  number 
of  joints,  the  technique  C *es  not  converge  to  a  unique  value  for  the 
joint  compliances.  The  table  below  shows  the  compliance  results 
obtained  considering  only  one  bending  mode  with  two  different  sets  of 
initial  compliance  estimates. 


It  is  seen  that  the  converged  solution  for  the  joint  compliance  is 
dependent  upon  the  initial  value  used  as  an  estimate  during  the  first 
iteration  when  the  number  of  modes  ueed  is  less  than  the  number  of 
unknown  joint  compliances. 

6.4  CURRENT  RESTRICTIONS 

The  joint  compliance  extraction  technique  currently  ha3  several 
restrictions.  One  of  the  restrictions  is  that  the  number  of  bending 
modes  used  must  equa.  or  exceed  the  number  of  unknown  joints  to  obtain 
accurate  results.  Also,  only  bending  cases  with  free-free  boundary 
conditions  can  be  run  using  the  joint  compliance  extraction  computer 
program,  and  a  method  of  handling  appendages  has  net  yet  been  devised. 

The  feasibility  of  the  joint  compliance  extraction  technique  has 
bean  demonstrated.  Several  limitations  exist  at  the  present  t.un&,  but 
it  is  felt  that  these  restrictions  can  be  overcome.  Note  has  been 
taken  of  a  premising  technique  for  handling  the  type  of  problem  being 
considered  here.  This  technique  has  recently  been  developed  by  Hall, 
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methc^  can  be  modified  and  applied  to  the  problem  of  e^acUna  loint 
compliances  from  measured  tactical  missile  elastic  mode  data.  §  J  * 
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Table  6-1 

SENSITIVITY  OF  JOINT  COMPLIANCE 
TO  INPUT  FREQUENCY  ERRORS 


Modal 

Frequency 

Variations 


;  C/C  Actual 

Using  Generated 

Modes 

Initial  Estimate 

.888 

.957 

.998 

.1 

1.06 

1.11 

.348 

.312 

.356 

.929 

1.26 

0 

1.24 

.866 

10. 

1.31 

1.33 

1.63 

.828 

1.13 

.946 

1.01 
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Table  6-2 


SENSITIVITY  OF  JOINT  COMPLIANCE 
TO  NUMBER  OF  MODES  USED 
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Section  /.0 

EXPERIMENTAL  VERIFICATION 


A  limited  test  program  using  highly  simplified  Joint  test  models 
has  been  conducted  to  illustrate  and  verify  seme  of  the  results  of  the 
parametric  analyses*  A  principal  source  of  airframe  Joint  compliance 
is  believed  to  result  from  strain  field  disturbances.  Almost  all 
Joints  by  definition  represent  a  significant  interruption  and  altera- 
tion  of  the  load  paths  between  sections.  In  mary  cases,  the  number  of 
load  paths  is  equivalent  to  the  number  of  principal  fasteners.  The 
conceptual  joint  model  discussed  in  Section  5.2.2  provides  a  useful 
design  basis  for  representing  both  axially  symmetric  and  non-axially 
symmetric  joint  test  specimens.  Physical  property  control  is  a  criti¬ 
cal  requirement  in  developing  simple  inexpensive  models  of  airframe 
joints.  The  design  approach  trken  and  described  in  this  section  has 
provided  a  total  of  12  test  specimens  in  three  control  groups,  having 
physical  properties  easily  defined  and  controlled,  and  requiring  only 
very  simple  machining. 

7.1  TEST  MODEL 

The  joint  simulation  models  consist  of  aluminum  tubing  with  circum¬ 
ferential  slots  cut  out  at  the  midspan  of  the  beams.  These  slots 
produce  local  strain  field  disturbances  under  loading  which  are  intended 
to  be  analogous  to  those  occurring  in  a  missile  airframe  joint.  Slot 
geometries  were  selected  to  provide  reductions  in  local  section  proper¬ 
ties  (El)  ranging  from  60  to  S7.5  percent  of  the  basic  cross  section. 

The  net  compliance  produced  by  these  slots  is  shown  to  be  significantly 
larger  t?.m  that  predicted  by  considering  only  the  area  reduction  over 
the  length  of  tha  slots.  The  strain  field  disturbance  is  propagated  a 
significant  distance  from  the  sectional  discontinuity.  In  deriving  an 
equivalent  beam  representation  for  the  test  models,  the  SI  distribution 
in  the  vicinity  of  the  joint  has  to  be  reduced  in  addition  to  the  El 
reduction  at  tha  joint,  or  a  local  joint  compliance  must  be  introduced 
at  the  joint  it;  order  to  predict  proper  mode  shape  and  frequency  data. 

The  joint  simulation  models  are  described  in  Figure  7-1.  Model 
number  1  is  a  uniform  tube  without  a  joint,  the  standard  or  reference 
for  this  set  of  models.  Models  2  thru  5  have  six  segments  of  material 
equally  spaced  around  the  circumference.  Model  numbers  6  thru  10  have 
three  equally  spaced  segments  of  material  remaining  from  the  slot  cuts. 

Model  numbers  11  and  12  are  designed  to  be  non-axially  symmetric  and 
produce  elastic  coupling  across  the  joint  by  the  mechanism  discussed  in 
Section  5.2.  The  segment  spacing  in  models  11  and  12  (8  *  40®,  20®) 
represents  an  estimated  elastic  coupling  of  45  and  85  percent  respectively. 
When  excited  Lu  either  of  the  two  planes  of  symmetry,  the  response  char¬ 
acteristics  will  show  only  one  peak  in  the  vicinity  of  the  first  mode.  When 
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tested  in  any  other  plane,  a  double  peaked  response  will  result.  The 
plane  with  maximum  cross  plane  coupling  will  result  at  45  degrees  from 
the  planes  of  symmetry. 

7.2  TEST  SETUP 

The  vibration  testing  of  the  joint  simulation  test  models  was  con¬ 
ducted  with  the  models  suspended  by  soft  springs  at  the  nodal  points  of 
the  first  bending  mode  as  shown  in  Figure  7-2.  This  figure  also  shows 
some  of  the  equipment  used  during  the  test.  A  miniature  accelerometer, 
shown  mounted  on  the  test  specimen,  and  a  Bentley  distance  detector 
were  used  to  make  response  measurements  on  the  test  models.  The  acceler¬ 
ometer  and  probe  outputs  were  fed  through  a  high  pass  filter  to  an 
oscilloscope  where  the  wave  form  was  monitored  and  to  a  counter  to 
measure  the  excitation  frequency.  An  oscillator  and  power  amplifier 
were  used  to  drive  an  induction  shaker  which  excited  the  test  models  in 
the  horizontal  plane. 

The  resonant  frequencies  of  the  joint  simulation  models  were  meas¬ 
ured  in  two  perpendicular  planes  for  the  symmetric  models  to  verify  that 
the  variations  in  frequency  due  to  manufacturing  and  measurement  toler¬ 
ances  were  small.  The  resonant  frequencies  of  model  number  11  and  12 
were  measured  in  the  two  planes  of  symnetry,  and  at  angles  of  thirty  and 
forty-five  degrees  frem  the  planes  of  symmetry. 

7.3  RESULTS 

The  results  of  the  vibration  testing  of  the  joint  simulation  models 
are  presented  in  Table  7-1.  With  the  exception  of  test  model  number  10, 
the  symmetric  joint  simulation  models  have  resonant  frequencies  which 
are  very  close  to  each  other  for  the  two  perpendicular  (A  and  B)  planes. 
With  specimen  number  10,  a  double  peaked  first  mode  was  observed. 
Obviously,  the  two  perpendicular  planes  tested  were  not  principal  stiff¬ 
ness  planes,  since  a  double  peaked  response  was  observed.  The  magnitude 
of  the  input  force  for  these  tests  was  so  small  that  response  amplitude 
measurements  could  net  be  made  off  resonance.  At  resonance,  very  little 
damping  existed,  and  the  peak  response  was  very  largo. 

Test  specimens  11  and  12  were  designed  to  produce  the  double  peaked 
first  mode  response  characteristics  discussed  in  Section  5.2.  Table  7-1 
shows  the  stiffness  ratios  in  the  two  principal  stiffness  planes  and  the 
resulting  resonant  frequencies.  The  first  mode  double  peaks  are  very 
widely  separated  (similar  to  the  type  of  response  shown  in  Figure  5rl6) 
due  to  the  differences  in  stiffness  between  the  two  planes. 

The  effective  compliance  produced  by  the  slot  discontinuities 
in  the  test  specimens  can  be  represented  by  either  of  the  joint  models 
described  in  Section  4.1.  Table  7r2A  lists  flexural  and  reduced  stiff - 
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ness  model  equivalences  for  each  test  specimen  based  on  matching  the 
measured  first  mode  frequencies.  In  the  first  case,  the  effective 
length  of  the  stiffness  discontinuity  induced  by  the  slot  is  determined, 
and  in  the  second  case,  the  rotational  compliance  of  the  equivalent 
joint  spring  is  determined.  Table  7**2B  shows  the  relationship  between 
test  specimen  effective  compliance  and  the  joint  "rating"  system  prev¬ 
iously  discussed  in  Section  4.2. 

A  plot  of  the  effective  length  of  the  stiffness  discontinuity 
expressed  in  test  specimen  body  diameters  is  shown  in  Figure  7-3.  As 
might  be  expected,  the  effective  length  decreases  as  the  number  of 
slots  (load  paths)  increase.  It  is  of  interest  to  note  that  for  the 
three  segment  models,  the  effective  length  ranges  from  about  0.2  tc  0.7 
diameters  with  the  average  of  approximately  0.5  falling  in  the  vicinity 
of  what  would  be  considered  a  good  to  moderate  joint. 

This  observation  is  perhaps  better  Illustrated  in  Figure  7-4, 
which  presents  test  specimen  geometric  stiffness  ratio  versus  frequency 
ratio.  The  numbers  adjacent  to  the  data  points  identify  the  test 
specimens.  The  6-segment  models  are  clearly  less  sensitive  than  the 
3 -segment  models  to  stiffness  discontinuities.  The  dashed  curve  in 
this  figure  represents  for  comparison  purposes  the  predicted  relation¬ 
ship  between  frequency  ratio  and  joint  stiffness  ratio  for  a  uniform 
beam  with  a  joint  at  the  mid-span.  This  curve  is  based  o.t  the  r.e3t 
specimen  fineness  ratio  (L/D)  of  15  and  a  reference  length  of  0.5  body- 
diameters  for  the  assumed  span  of  the  local  stiffness  reduction.  Test 
points  close  to  the  dashed  curve  imply  effective  lengths  in  stiffness 
loss  close  to  the  0.5  diameter  assumption.  The  branch  curve  passing 
through  point  8  connecting  data  for  test  specimens  11  and  12  reveals 
the  dramatic  influence  of  unsymmetric  segment  spacing.  In  this  instance, 
the  segment  areas  for  the  three  specimens  (8,  11,  12)  are  the  same  and 
it  is  interesting  to  note  that  toe  variation  with  segment  spacing  approaches 
the  dashed  curve  quite  closely. 

7.4  CORRELATION  ANALYSIS 

In  an  effort  to  obtain  a  relatively  accurate  analytic  determination 
of  the  equivalent  beam  stiffness  for  comparison  with  test  results  a 
relatively  detailed  finite  element  analysis  was  undertaken  for  one  test 
specimen  -  number  10.  This  analysis  was  performed  using  the*  STARDYNE 
computer  program,  Reference  8. 

The  beam  with  a  concentric  circular  cross  section  is  interpreted 
as  a  cylindrical  shell  with  a  length  equal  to  four  diameters.  The  notch¬ 
ing  is  placed  at  the  mid -section.  The  ends  of  the  shell  are  given  a 
stress  distr ’.Dution  which  varies  linearly  in  one  direction  and  which 
results  in  a  uniform  moment  distribution  along  the  axis  of  the  associated 
beam.  The  shell  and  loading  are  shewn  in  Figure  7-5a. 
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Due  to  symmetry  only  one  fourth  of  the  geometric  configuration 
need  he  modeled  as  shown  in  Figure  ?-5b.  The  finite  element  analysts 
used  described  this  portion  of  the  geometric  configuration  with  a  set 
of  triangular  finite  elements.  This  element  accounts  for  both  shell 
msmbi  -e  and  bending  stresses.  The  distribution  of  finite  elements 
used  is  shown  in  Figure  7-6.  This  figure  depicts  the  element  distribu¬ 
tion  for  the  unfolded  shell  segment  so  that  it  can  be  placed  on  a  plane. 

The  model  used  has  480  triangular  elements  ami  275  node  points. 
Since  there  are  6  degrees  of  freedom  for  each  node  point  the  model  has 
1650  degrees  of  freedom.  The  application  of  boundary  conditions  for 
the  quarter  segment  requires  the  constraining  of  108  degrees  of  freedom, 
lhus  the  model  has  1542  unconstrained  degrees  of  freedom. 


fiie  analysis  yielded  a  set  of  node  displacements  and  rotations. 
This  in  turn  was  used  to  determine  an  equivalent  compliance  value  for 
spring  placed  at  the  notched  section.  The  value  of  the  computed 
compliance  is  7.03  :<  (10)"6  rad/in-lb.  The  measured  frequency  for 
model  number  10  Implies  a  compliance  value  of  7.6(10)“6  rad/in-lb.  The 
finite  element  analysis  compliance  value  agrees  with  the  eyperim°ntal 
■alue  within  7.5  per  cent. 
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Table  7-1 


Vibration  Test  Results 
Free -Free  Joint  Simulation  Models 


Specimen 

Number 

■3 

** 

(Hz) 

V* 

r— - 

V* 

|Qgj| 

1 

526. 

526. 

1.0 

1.0 

1.0 

2 

511. 

510. 

.971 

.970 

0.4 

3 

492. 

492. 

.935 

.935 

0.2 

4 

476. 

475. 

.905 

.903 

0.1 

5 

453. 

454. 

.861 

.863 

0.05 

6 

484. 

484. 

.920 

.920 

0.4 

7 

436. 

436. 

.829 

.829 

0.2 

8 

403. 

402. 

.766 

.764 

0.1 

9 

382. 

380. 

.726 

.722 

0.05 

10 

356. 

356. 

.677 

.677 

0.025 

360. 

360. 

.684 

.684 

MB  •« 

11 

334. 

451. 

.635 

.857 

.0554/. 1446 

12 

236. 

472. 

.449 

.897 

.0162/. 184 

■£  «  First  bending  mode  frequency  of  reference  specimen  (No.  1). 

£  s  First  bending  mode  frequency  in  the  A  plane. 

$6  -  First  bending  mode  frequency  in  the  B  plane. 

=  Sectional  stiffness  of  slotted  region  of  the  tuba. 

-  Sectional  stiffness  of  the  basic  tube. 
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Table  7-24 

COMPLIANCE  AND  COMPLIANCE  RATIO  RESULTS 
FOP  JOINT  SIMULATION  MODELS 


— 1 —  ■  ■ — ' - 

Specimen 

No.  of 

Freq. 

Geometric 

Effective 

Effective 

Number 

Seg. 

Ratio 

Stiffness 

Length 

Compliance 

Ratio,  IR 

//D 

Rad/ln-LB 

(10)6 

1 

1.0 

1.0 

0 

2 

6 

.571 

0.40 

.245 

.4 

3 

G 

.935 

0.20 

.230 

1.0 

4 

6 

,905 

0.10 

.150 

1.4 

5 

i 

6 

.861 

.05 

.106 

2.2 

6 

3 

.520 

0.40 

734 

1.2 

7 

3 

.829 

0.20 

.666 

2.9 

8 

3 

.766 

0.10 

.459 

4.5 

3 

.726 

0.05 

.274 

5.7 

3 

.677 

0.025 

.179 

7.6 

3 

.857 

.145 

.359 

2.3 

3 

.635 

.055 

.516 

9.6 

3 

.897 

.184 

.332 

1.6 

3 

.449 

.016 

.336 

25.5 

Table  7-2B 

TEST  SPECIMEN  EFFECTIVE  COMPLIANCE  COMPARISON 
WITH  JOINT  RATING  CLASSIFICATION 


mm 

Joint  "Rating" 

Excellent 

Good 

Moderate 

Ce(10)6  Rad/ln-Lb 

0.1 

1.0 

1G.0 

flft 
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Section  8.0 
CONCLUSIONS 


In  reviewing  the  results  of  this  initial  study  of  the  structural 
dynamic  properties  of  tactical  missile  joints  the  major  points  would 
appear  to  be  the  following. 

(1)  The  most  significant  directly  related  work  to  date  is  that  of 
Alley  and  Leadbetter  (Reference  1)  in  which  the  compliance 
characteristics  of  various  joints  are  classified  from  "excel lent*' 
(small  compliance)  to  "loose". 

(2)  The  compliance  and  damping  characteristics  of  typical  tactical 
missile  joints  can  have  a  powerful  influence  on  airframe 
structural  dynamic  response.  A  single  joint  of  "moderate" 
compliance  placed  at  the  midspan  of  a  typical  missile;  for 
example,  could  be  expected  to  reduce  the  first  mode  frequency 
by  approximately  35  percent. 

(3)  The  structural  efficiency  of  typical  missile  joints  judged  on 
the  basis  of  maintaining  airframe  stiffness  through  the  local 
region  of  the  joint  is  extraordinarily  low  with  a  "good"  joint 
rated  near  40  percent,  a  "moderate"  joint  near  5  percent,  and 
a  loose  joint  generally  below  1  percent  (efficiency  in  this 
sense  equals  percent  of  local  stiffness  retained  over  a  length 
of  one-half  a  body  diameter). 

(4)  Results  obtained  with  simple  test  models  suggest  that  the 
number  of  load  paths  and  their  spacing  around  the  periphery  of 
the  joint  are  the  determining  parameters  in  joint  compliance. 

The  poor  performance  of  some  types  of  ring  joints  is  suspected 
to  be  traceable  to  relatively  low  axial  preloads  (although 
assembly  torques  may  be  high  due  to  friction)  and  a  strong 
likelihood  of  only  three  points  of  contact  unequally  spaced. 

(5)  The  role  of  unequal  load  path  spacing  in  joints  as  a  source 
of  elastic  coupling  has  been  illustrated  both  by  analysis  and 
test,  and  a  simple  conceptual  model  presented  which  is  believed 
to  offer  a  plausible  explanation  and  insight  for  this  mechanism. 

(6)  A  technique  for  extracting  the  joint  compliance  values  from  a 
set  of  measured  missile  elastic  mode  frequencies  and  shapes 
has  been  developed  to  a  preliminary  stage.  The  technique, 
which  is  based  on  the  method  of  steepest  ascent  used  originally 
in  trajectory  optimisation,  has  been  shown  to  be  feasible  for 
the  present  application.  Certain  limitations  are  present  in 
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the  h..  nod  as  currently  implemented.  These  limitations  ate 
not  inherent  in  the  technique.  The  removal  of  these  limitations 
will  result  in  a  powerful  and  useful  analytical  tool  for  the 
determination  of  flexural  joint  compliances  from  measured 
missile  elastic  mode  test  data. 

(7)  Finite  element  structural  analysis  approaches  are  believed  to 
offer  excellent  potential  as  analytical  tools  in  estimating 
joint  load/deflection  characteristics.  Selection  of  the  model 
and  the  interface  load  path  assumption  are  critical 
considerations . 

The  next  phase  of  this  study,  identified  as  "Pha9e  2,  Data 
Acquisition  and  Analysis  Extension  Phase"  is  directed  at  accomplishing 
the  following  objectives. 

(1)  Developing  the  analytical  techniques,  which  were  formulated 
in  Phase  1,  into  operational  methods. 

(2)  Applying  these  operational  methods  to  the  data  collected  in 
the  Phase  1  study. 

(3)  Developing  approaches  to  modeling  the  more  significant  damping, 
cross  coupling  and  nonlinear  stiffness  joLnt  characte  )tics. 

(A)  Conducting  a  limited  test  on  full  scale  missile  joint  hardware. 

(5)  Performing  a  correlation  analysis  to  provide  a  basis  for 

verification  of  the  analytical  techniques  developed  in  Phase 

2. 
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APPENDIX  T 

AEROSPACE  INDUSTRY  SURVEY 


An  aerospace  Industry  purvey  on  structural  dynamic  properties  of 
missile  airframe  joints  has  been  conducted  in  the  form  of  a  questionnaire. 
This  questionnaire  was  distributed  with  the  Minutes  of  the  Aerospace 
Flutter  and  Dynamics  Council  Meeting  held  May  14  -  15,  1969  in  San  Antonio, 
Texas.  The  intent  >t  the  survey  was  to  gather  information  from  which  a 
list  of  joints  in  common  usage  and  the  structural  dynamic  characteristics 
of  importance  for  joints  could  be  compiled.  The  primary  results  of  the 
aerospace  industry  survey  are  presented  in  the  form  of  a  table  (Table  £-1). 


The  authors  wish  to  express  their  appreciation  to  the  following 
individuals  and  organizations  woo  have  provided  the  information  given  in 


this  Appendix. 

i'r.  Thomas  W.  Miller 

The  Hoeing  Company 

Seattle,  Washington 

Mr.  I.  baker 

Mr.  F.  A.  Figge 

Mr.  R.  J.  Oedy 

Hughes  Aircraft  Co. 

Canoga  Park,  California 

Mr.  Bet «.  M.  Hall 

Mr.  Gerry  Kahre 

McDonnell  Ocuglas  Astronautics  Co. 

Santa  Monica-  California 

Mr.  Craig  S.  Porter 

Mr.  William  J.  Werback 

Naval  Weapons  Center 

China  Lake,  California 

Mr.  H.  M.  Marshall 

Philco  Ford  Corp. 

Newport  Beach,  California 

Table  l-l  shows  specific  missile  and  joint  properties  obtained  from 
the  aerospace  Industry  survey  plus  data  obtained  here  at  the  Pomona 
Operation  of  General  dynamics.  Only  a  limited  number  of  joint  compliance 
values  have  been  obtained  from  the  industry  survey,  and  all  of  the  com¬ 
pliances  have  been  flexural  compliances.  No  shear  compliances  have  been 
reported,  supporting  the  premise  that  flexural  compliance  is  dominant  to 
shear  compliance  in  importance  to  structural  dynamic  analyses.  Not  all 
of  those  responding  to  the  questionnaire  represented  joints  as  discrete 
springs,  but  rather  seme  considered  the  bending  stiffness  to  be  some 
reduced  value  of  stiffness  in  the  area  of  the  joint.  The  effect  of  the 
airframe  joints  on  the  first  bending  mode  frequency  has  been  obtained  for 

six  tactical  missiles.  The  results  are  presented  in  Table  1-2. 

A  decrease  in  airframe  stiffness  as  the  number  of  joints  is  increased  can 
be  seen. 
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APPENDIX  II 
3IBLI0GRAPHY 


A  literature  search  was  conducted  on  the  topic  of  the  structural 
dynamic  properties  of  tactical  missile  airframe  joints.  The  following 
literature  sources  were  used. 

1.  Engineering  Index,  1960  to  August  1969 

2.  Scientific  and  Technical  Aerospace  Reports,  1960  to  1969 

3.  Defense  Documentation  Center  search  using  the  following 
search  words;  missile  joints,  joints,  mechanical  joints, 
joint  stiffness,  joint  flexibility,  and  dynamic  properties 
of  joints. 

A  list  of  publications  which  contain  information  relative  to  the 
topic  follows  with  a  brief  summary  of  each  of  the  articles. 

II-l  Alley  Jr.,  V.  L.  an<*  ’.eadbetter,  S.  A.,  "The  Prediction  and 
Measurement  of  Natural  Vibrations  on  Multistage  Launch  Vehicles*', 

American  Rocket  Society  Launch  Vehicles:  Structures  and  Materials 
Conference  Report  April  1962. 

Results  of  an  analytical  and  experimental  study  to  determine 
natural  frequencies  of  a  multistage  research  rocket  are  presented. 

Also  studies  of  the  effects  of  mechanical  joint  looseness  and  an  empir¬ 
ical  treatment  or.  joint  flexibility  are  discussed.  Some  typical  missile 
joints  tre  illustrated  and  classified  from  excellent  to  loose  based  on 
joint  compliance  and  diameter. 

II-2  Barton,  M.  V.,  "important  Research  Problems  in  Missile  and  Space¬ 
craft  Structural  Dynamics",  NASA  :.’N  D-X296,  May  1962. 

The  author  provides  a  broad  overview  of  structural  dynamic  problems 
of  missile  and  space  vehicles  and  attempts  to  identify  those  areas  ir. 
which  further  research  would  be  fruitful  in  enhancing  the  state-of-the- 
art.  Although  this  report  was  written  nearly  10  years  ago  and  much  pro* 
gress  has  been  made  since  then,  much  of  what  is  indicated  in  this  report 
is  pertinent  today.  He  does  not  explicitly  identify  missile  joints  as 
a  problem  area  although  this  is  implicit  in  several  of  the  problem  areas 
identified. 

XI-3  Collins,  J.  D.,  and  Thomson,  W.  T.,  "The  Eigenvalue  Problem  for 
Structural  Systems  with  Statistical  Properties",  AIAA  Journal,  Vol.  7, 

No.  4,  April  1969,  pp  642-648. 
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The  treatment  of  matrix  eigenvalue  problems  for  structural  sys- 
terns  with  statistical  properties  Is  presented.  Using  the  techniques 
described  in  the  paper,  Eiger.vrlues  can  be  determined  in  a  statistical 
form  for  a  missile  which  has  its  joint  stiffness  properties  described 
in  a  statistical  form.  Variations  which  can  occur  in  serial  production 
can  be  accounted  for  and  consequently  the  technique  could  prove  to  be 
very  useful. 

II-4  DeVries,  G.,  "Design  of  Joints  for  Cylindrical  Sections",  U.  S. 
Naval  Ordnance  Test  Station,  China  Lake,  California,  Report  No.  IDP807, 
April  1960. 

The  author  classifies  the  various  types  of  circular  joints  used 
in  missile  applications  and  presents  a  set  of  criterion  for  evaluating 
the  desirability  of  the  various  classes.  The  criteria  used  do  not 
include  structural  dynamic  considerations. 

II-5  Grarabell,  II.  V.,  "A  Compendium  of  Structural  Joints  for  Assembly, 
Field  and  Flight  Separation  on  Missiles*,  Boeing  Co.,  Seattle,  Report 
No.  D2 -12591 1-1-Rev-l.cr-B,  July  68. 

The  author  assembles  a  cross-section  of  state-of-the-art  designs 
of  missile  joints  used  by  the  3oeing  Company.  The  merits  of  the  various 
joints  are  evaluated,  however,  structural  dye  tmic  response  considera¬ 
tions  are  not  examined. 

II-6  Hanks,  B.  R.  and  Stephens,  D.  G.,  Mechanisms  and  Scaling  of 
Damping  in  a  Practical  Structural  Joint”,  Shock  and  Vibration  Bulletin 
No.  36  Part  4,  pp  1  -  8,  January  1967. 

The  authors  report  on  an  investigation  directed  at  determining 
the  effect  of  geometric  scale  on  the  damping  of  a  beam  joint  assembly. 

In  essence  the  decay  of  the  fundamental  mode  of  four  similar  cantilever 
configurations,  varying  in  scale  from  20  to  l,  were  experimentally 
determined  for  various  bolt  tension  loads.  It  is  concluded  that  damping 
is  inversely  proportion'll  to  model  size. 

11-7  Kalinia,  N.  G.,  Lebedeve,  Y.  A.,  ei  al.,  "Structural  Damping  in 
Permanent  Joints",  Translation  Division  WP-AFB  Ohio,  FID-TT-63-755/1  +  2, 
May  64. 


The  authors  report  on  various  investigations  on  structural 
damping  and  consequently  the  studies  are  directed  at  damping  mechanisms 
associated  with  structural  joints.  The  authors  identify  the  dominant 
damping  mechanism  associated  with  bolted  and  riveted  joints  as  coulomb 
damping  associated  with  interface  frictional  forces.  It.  should  be 
pointed  out  that  more  recent  studies  put  forward  evidence  indicating 
that  this  type  of  damping  is  in  fact  viscous  damping  associated  with  the 
talent iai.  movement  of  air  at  the  interface. 
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T.I-8  Maidanik ,  G.  and  Ungar,  E,  E.,  “Panel  Loss  Factors  Due  to  Gas- 
Pumping  at  Structural  Joints'1,  NASA  CR-954,  November  1967. 

The  authors  report  on  a  theoretical  and  experimental  study  cf.  gas 
pumping  in  riveted  joints.  Results  from  a  mathematical  model  are 
compared  with  experimental  results.  The  results  indicate  that  for  the 
higher  frequencies  gas  pumping  is  the  dominant  mechanism.  The  authors 
also  note  that  structural  damping  measurements  conducted  at  an  ambient 
ores sure  may  yield  misleading  results  for  structures  which  operate  in 
a  rarefied  atmosphere. 

II-9  Mead,  D.  J.,  and  Eaton,  D,  C.  G.,  “Interface  Damping  at  Riveted 
Joints,  Part  I  -  Theoretical  Analysis",  Wright -Patterson  Air  Force  Base, 
Ohio,  ASD  Technical  Report  61-467  Sept.  1961. 

A.  theo .  etical  examination  is  made  of  damping  in  a  riveted  lap 
joint  having  a  viscc-elastic  interfacial  layer.  Simple  design  rules 
are  given  for  maximum  interface  damping. 

11-10  Mentel,  T.  L. *  “Joint  Interface  Layer  Damping,"  Transactions  of 
ASME  Vol.  89,  Series  B,  No.  4,  Nov.  1967. 

An  experimental  and  theoretical  study  of  damping  in  thin  visco¬ 
elastic  layers  is  presented.  The  author  notes  that  the  damping  mechanism 
is  dominated  by  the  analytically  simple  shear  mechanism. 

II-ll  Rubenstein,  N.,  Sigillito,  V,  G.,  and  Stadter,  J.  T.,  “Upper  and 
Lower  Bounds  to  Bending  Frequencies  of  Non-uniform  Shafts,  and  Applica¬ 
tions  to  Missiles,"  The  Shock  and  Vibration  Bulletin  No.  38,  Part  2, 

Aug.  1968,  pp  169  -  176. 

The  authors  present  a  method  for  computing  upper  and  lower  bounds 
to  bending  frequencies  of  non-uniform  shafts  which  have  their  El  distri¬ 
bution  precisely  defined.  The  authors  indicate  that  this  type  of  determi¬ 
nation  is  important  in  missile  developments  and  undoubtedly  it  is, 
however,  precise  definition  of  El  distributions  for  missiles  is  difficult 
to  establish  due  primarily  to  joint  behavior.  Consequently  the  bourds 
are  not  necessarily  meaningful  for  many  missile  applications. 

11-12  Smith,  F.  A.,  "Acoustic  Response  Analysis  of  Large  Structures," 

The  Shock  and  Vibration  Bulletin,  No.  39,  Part  3,  Jan.  1969,  pp  55  -  64, 

The  technical  approach  and  results  of  an  acoustic  response  analysis 
up  to  200  Rz  for  a  large  complex  structure  is  presented.  An  elastic 
structural  model  which  used  3000  degrees  of  freedom  and  based  on  beam  and 
plate  elements  together  with  some  one  dimensional  elements  was  used.  A 
surprisingly  good  correlation  with  measured  results  was  obtained.  This 
type  of  detailed  analysis  is  desirable  because  cf  the  quality  of  the 
results  which  can  be  obtained.  An  analysis  with  this  degree  of  complex - 
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Ity  (for  r<my  missile  applications)  is  often  limited  due  to  an  inability 
to  describe  joint  behavior. 

11-13  Trotter,  W.  D.,  Rauch,  G.  G.  and  Muth,  B.  V...  'Tlissile  Dynamic 
Response  Tests,*'  Boeing  Company,  Seattle,  Report  No.  T2AGM20308-12,  l, 
-3-69. 


Results  for  the  SRAM  missile  structural  dynamic  response  tests 
are  presented.  The  effect  of  the  compliance  of  the  airframe  joints  on 
the  measured  bending  mode  shapes  is  readily  apparent. 

11-14  Ungar,  E.  E.  and  Carbonell ,  J.  R. ,  "On  Panel  Vibration  Damping 
Due  to  Structural  Joints,"  A1AA  Journal,  August  1966,  pp  1385  -  1390. 

The  author  reports  on  an  experimental  study  of  damping  mechanisms 
for  bolted  joints.  He  concludes  that  high  frequency  damping  is  primarily 
due  to  the  pumping  of  air  resulting  from  surfaces  moving  away  from  and 
toward  each  other. 
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